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Abstract. We study theoretically the dynamics of living polymers which can add and subtract monomer 
units at their live chain ends. The classic example is ionic living polymerization. In equilibrium, a delicate 
balance is maintained in which each initiated chain has a very small negative average growth rate ( "veloc- 
ity") just sufficient to negate the effect of growth rate fluctuations. This leads to an exponential molecular 
weight distribution (MWD) with mean N. After a small perturbation of relative amplitude e, e.g. a small 
temperature jump, this balance is destroyed: the velocity acquires a boost greatly exceeding its tiny equi- 
librium value. For e > ec ~ l/N^^^ the response has 3 stages: (1) Coherent chain growth or shrinkage, 
leaving a highly non-linear hole or peak in the MWD at small chain lengths. During this episode, lasting 
time Tfast ~ .^V, the MWD's first moment and monomer concentration m relax very close to equilibrium. 
(2) Hole-filling (or peak decay) after tah ~ e^N^. The absence or surfeit of small chains is erased. (3) 
Global MWD shape relaxation after TsIow ~ -^V^. By this time second and higher MWD moments have 
relaxed. During episodes (2) and (3) the fast variables [N , m) are enslaved to the slowly varying number 
of free initiators (chains of zero length) . Thus fast variables are quasi-statically fine-tuned to equilibrium. 
The outstanding feature of these dynamics is their ultrasensitivity: despite the perturbation's linearity, the 
response is non-linear until the late episode (3). For very small perturbations, e < Ec, response remains 
non- linear but with a less dramatic peak or hole during episode (1). Our predictions are in agreement with 
viscosity measurements on the most widely studied system, Q-methylstyrene. 

PACS. 82.35.-X Polymers: properties; reactions; polymerization. 05. 40. -a Fluctuation phenomena, random 
processes, noise, and Brownian Motion. 87.15.Rn Biomolecules: structure and physical properties; Reactions 
and kinetics; polymerization. 
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Living polymers are intriguing examples of soft matter spond dynamically via polymerization and depolymeriza- 

and have major roles in technology and biology. The clas- tion reactions allowing a new thermodynamic equilibrium 

sic system, of major importance for the synthesis of high- to be attained. These systems are "living" in the sense 

grade polymer materials and the subject of a long his- that a change in their environment leads to a new equi- 

tory of fundamental experimental and theoretical study, librium molecular weight distribution (MWD). All of this 

is living anionic polymerization jHEllS ,4,5,6, 111111111101 should be compared to conventional inert polymer mate- 

|TT|[T^|T^[ni|T^P^|T7|[TBl[T^|^|niP2H23„24t|^ This rials whose MWDs are frozen, 
class of polymerization is widely employed to manufac- 
ture polymers with nearly monodisperse MWDs and con- 

trolled architectures such as block and star copolymers The subject of this paper is the dynamical response 
Eg. A recent variant on this theme with immense poten- of ^'""'"^S polymers. We study the class of hvmg polymers 
tial technological impact is living free radical polymeriza- ^^^^^ exemplified by hvmg ionic polymerization, having 
tion 27,28 . Other examples include worm-like surfactant following characteristics: (i) monomers add and sub- 
micelles |2Hl3()l;-!ll;-i2l33l:^4l:ihl:ffi| and biological polymers ^ract from chain ends and (u) the total number of chains 
such as actin and tubulin filaments |S71EHlE3EnilIIllll ^^'^ biological world, actm and microtubule 
whose special properties are exploited by living cells for filaments satisfy condition (i), and the dynamics of both 
motility and structural integrity filaments and filament caps are frequently described by 

As polymeric or pofymer-like materials, the unique situation (ii) m Wormlike micelles on the other hand 
feature of these systems is that the chains are dynamic ^O":!*; ^^^ses grow from their ends only |ini|3D| (condi- 
tion (i)) but do not satisfy condition (ii): the number of 

e-mail: bo8@columbia.edu micelles is not conserved. Interestingly, this makes their 

^ e-mail: dv35@columbia.edu dynamical response very different (see discussion section). 
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Fig. 1. (a) Schematic of initiation reaction in ionic polymer- 
ization. Empty (filled) circle represents initiator (monomer), 
(b) Polymerization and depolymerization reactions at the ionic 
chain end of a polymer chain of A*' units. 

Our presentation will focus on living ionic polymer- 
ization, though most of our results are completely gen- 
eral. In these systems, polymers grovif reversibly from their 
ionic ends (see fig. Due to their "living" nature, if 
fresh monomer is added polymers will grow until the extra 
monomer has been consumed; further addition of a dif- 
ferent monomer species, say, will resume polymer growth. 
The main achievement of this technique is a distribution of 
chain lengths with very small Poisson fluctuations around 
the mean. This monodisperse MWD is not however a true 
equilibrium distribution. For sufficiently long times slow 
depolymerization reactions become important and lead to 
a true equilibrium MWD which is in fact extremely broad. 
For many applications, for example styrene, the timescale 
to reach equilibrium is extremely large and depolymeri- 
zation-induced broadening effects are unimportant 4A. 
However, in some cases equilibration times are accessi- 
ble, an example being a-methylstyrene whose equilibrium 
properties have been studied in many experiments 
1^,6,7,8,9,10,11,12,13,. 

The reason that chain length fiuctuations in equilib- 
rium are so large is that even when a state is reached where 
the monomer concentration reaches a steady state value 
such that polymer growth rates exactly balance depoly- 
merization, through polymerization and depolymerization 
reactions monomers continue to be reshuffled among liv- 
ing chains. Thus even if one starts with an essentially 
monodisperse MWD as shown in fig. El^a), random ad- 
dition and subtraction of monomers from chain ends will 
result in a "diffusive" random walk motion of chain ends in 
A'^-space (where N is the number of monomer units added 
to initiator) with a certain characteristic "diffusivity" D. 
Since there is no restoring force to this motion, chain ends 
will eventually diffuse to distances of the same order as 
the mean chain length, resulting in a broad distribution 
as shown in fig.EIb). In a simple mean-field picture where 
monomer-monomer excluded volume interactions are ne- 
glected this leads to the Flory-Schultz equilibrium distri- 
bution gniiiniEi 

0oo(iV) = ^e-^/^- , (1) 

where N^o is the equilibrium mean length. 

The thermodynamic properties of equilibrium living 
polymers have been studied in a large number of experi- 
mental and theoretical works. According to theory, going 




Fig. 2. Fluctuations of living polymer chain length in equilib- 
rium are very large. For clarity, living chains are depicted as if 
extending from an imaginary substrate. Living ends are shown 
as filled circles. Even if one starts with a nearly monodisperse 
MWD as in (a) and with monomer concentration (grey back- 
ground) such that polymerization and depolymerization rates 
balance, random polymerization and depolymerization reac- 
tions eventually broaden the distribution up to distances of 
the same order as the mean A'^, as shown in (b). Note the re- 
sulting MWD includes chains of "zero" length (free initiators). 



beyond mean field by accounting for excluded- volume in- 
teractions leads to power law modifications jl8lll9lE(11l34| 
to eq. Moreover, an analogy jl4lll5l[TI)lll7| between a 
living polymer system and the much studied Ising spin sys- 
tem in the formal limit of spin dimensionality n becoming 
zero (similar to the well-known analogy in the thermody- 
namic theory of semidilute polymer solutions |47l48j which 
can in fact be formulated equivalently to a living polymer 
system showed that the polymerization phase transi- 
tion near the polymerization temperature is second-order 
with a number of characteristic power-law exponents. The 
results of a large number of experiments by the group of 
Greer, measuring both the MWD itself 13 and other ther- 
modynamic properties 2 4 are apparently consistent with 
the above theories, favoring somewhat the non mean-field 
approach. 

The most essential feature of living polymers however, 
namely their ability to respond dynamically to external 
perturbations, remains largely unexplored. How does an 
equilibrium living polymer system respond to a small per- 
turbation? A conceptually simple perturbation is the addi- 
tion of a small amount of extra monomer (an "m-jump"). 
A similar perturbation, easier to perform in practice, is a 
"T-jump," i. e. a small sudden change in temperature. Un- 
der such a perturbation, how rapidly and in what manner 
will the monomer concentration, MWD and mean chain 
length reach a new equilibrium? In this paper we address 
these questions theoretically. We will compare our pre- 
dictions in section 8 with small T-jump viscosity relax- 
ation experiments by Greer et al. |11) . This same group 
has also monitored monomer and MWD dynamics fol- 
lowing much stronger perturbations, finite temperature 
quenches [T^IOj (see fig.OJ. 

Our main broad conclusion is that, counter- intuitively, 
the response to a very small perturbation is extremely 
strong: even though the initial and final MWDs (the latter 
attained after relaxation is complete) are very close to one 
another, the MWD at intermediate times is very strongly 
perturbed from equilibrium, i. e. its shape is very different 
from either the initial or final equilibrium shapes. In this 
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Fig. 3. Schematic of one of the classes of experiment performed 
by Greer's group [T^IT^ . The mix of monomer, initiator and 
solvent is stored for sufficiently long above the polymerization 
transition temperature Tp such that all initiator is activated. At 
this point all living chains are of "zero" length; we refer to these 
as "free initiators" (bifunctional dimers in this case). After 
a quench below Tp, polymerization onsets and the system's 
response to abrupt temperature changes is then studied. In 
this paper we study more elementary perturbations: very small 
T-jumps of an already equilibrated system. 



sense, there are no small perturbations: a small thermo- 
dynamic perturbation well-described by a linear suscepti- 
bility (i.e. one inducing a small change in the equilibrium 
state) has nonetheless a large dynamical effect; intermedi- 
ate states deviate strongly from equilibrium in the sense 
that observables are perturbed in a highly non-linear man- 
ner. We refer to this as dynamical ultrasensitivity. 

A hint of this sensitivity is already apparent in equilib- 
rium, where a very delicate balance is established between 
growth and shrinkage such that the net polymerization ve- 
locity V is essentially zero, of order 1/ N^o (due to the small 
fraction of special non-depolymerizable chains of length 
TV = 0). Now the effect of a small T-jump (say) is to pro- 
duce a small but finite velocity which overwhelms this tiny 
velocity value and destroys the delicate balance sustain- 
ing equilibrium. Thus if for example v becomes positive, 
the MWD will uniformly translate toward larger molecular 
weights, creating a depletion in the MWD at small chain 
lengths, as shown schematically in fig0| This translation 
process stops only after sufficient monomer has been poly- 
merized such that the monomer concentration drops close 
to its destination equilibrium value, corresponding to a re- 
covery of the zero net growth rate equilibrium condition. 
We find that only on a much longer timescale does "dif- 
fusion" of chain ends fill up the hole and lead the system 
to the destination equilibrium MWD, as shown in fig. 01 

Why does this hole develop in the MWD? The rea- 
son is that chain length scales homogenized by diffusion, 
(-Dt)"'^/^, are smaller than scales affected by translational 
motion, vt, for times longer than the timescale t* « Djv^ . 
Thus if t* is shorter than the timescale needed for transla- 
tion to stop, a depletion hole develops in the MWD. But 
even when translation is completed earlier than t*, we 
find that the perturbation is still nonlinear; even though 
no hole is formed there is still a large amplitude reduction 
in the MWD, much larger than the relative magnitude of 
the initial perturbation, e. 



This work is related to earlier theory of Miyake and 
Stockmayer |44j who, following a prior treatment by Brown 
and Szwarc |50| . studied analytically and numerically the 
dynamics in the special case where initially all monomer 
is unpolymerized (see fig. [SJ. Their analytical results, ap- 
plicable in the case of very small depolymerization rates, 
showed a 2-stage relaxation. During the first, lasting a 
time ~ iVoo , living chains grow coherently from zero length 
up to the equilibrium value. In the second stage the shape 
of the sharply peaked MWD relaxes slowly toward the 
equilibrium exponential distribution. Although their per- 
turbation analysis broke down during this stage, Miyake 
and Stockmayer estimated that its duration scales as . 
Later theoretical studies addressed the second stage dy- 
namics 51,52, 53 of this special case. 

A short announcement of the present work has al- 
ready appeared jSl]. The structure of the paper is as 
follows. In section 2 we establish the dynamical equa- 
tions obeyed by living polymers which show the nonlinear 
coupling between the monomer population and the living 
chain MWD. For simplicity, we neglect excluded volume 
interactions and chain length dependent polymerization 
rates. Starting from the equations of section 2 we then 
analyze the response to a T-jump. This response depends 
on the sign of the velocity induced by the perturbation. 
The case of positive initial velocity is analyzed in sections 
3 and 4, while the negative in section 5. The special case 
of very small T-jumps is analyzed in section 6. In section 
7 we show that the results of the previous sections di- 
rectly generalize to small perturbations of arbitrary form. 
Finally we conclude with a discussion of the results, and 
the experimental outlook. 

2 Living Polymer Dynamics: 
Monomer-Polymer Coupling 

An important aspect of living polymers dynamics is the 
coupling between the two distinct species present in the 
solution, free monomers and living chain MWD: on the 
one hand living chain growth rates depend on the concen- 
tration of free monomer and on the other hand, due to 
mass conservation, monomer concentration if a function 
of the living chain MWD. Stated differently, living chains 
grow according to a velocity field which is self-consistently 
updated depending on the response of the MWD. In this 
section we develop the equations obeyed by monomer and 
MWD. The dynamical response to perturbations is dis- 
cussed in the following sections. 

A crucial feature of the living polymer systems we 
study is that the number of living chains remains fixed. 
This is what arises in the experiments of Greer et al. 
OSEli where this number is determined by the amount 
of initiator initially added in the solution. Denoting zq the 
concentration of living chains (i. e. chains which have been 
initiated), mass conservation implies that monomer con- 
centration, mj, and MWD, (f>t{N), [t denotes time) obey 

/•oo 

mt = mtot - i^Nt , Nt= dN N (pt{N) , (2) 

"'0 



Ben O'Shaughnessy, Dimitrios Vavylonis: Dynamics of Living Polymers 





eN 



Fig. 4. There are three stages in the relaxation of living polymers to equilibrium. During the first, living polymers grow 
coherently as in (a). The translation stops when they have consumed enough monomer such that monomer concentration 
drops to a value for which polymerization and depolymerization balance (shown in (b)). During the second stage, (c), through 
"diffusive" random polymerization and depolymerization reactions, the MWD slowly fills the depletion of chains of short lengths. 
The third and final stage entails global shape relaxation of the MWD on a difi'usion timescale corresponding to Noo ■ 



where mtot is the total monomer concentration (includ- 
ing polymerized monomers), and Nt is mean chain length. 
Here (j)t is normalized to unity and N — corresponds to 
free initiator. The coupling between monomer concentra- 
tion and the MWD is manifest in eq. 

Now the dynamical equations equations obeyed by (j)t 

are 



dMN) _ djt{N) 



dt 
Jt{N) 



dN 

VtMN) - 



ON 



Jt(0) = 0, (3) 



where 



rrit 



{k+mt + v-)/2 . (4) 



Here fc+ is the propagation rate constant and is the 
depolymerization rate. The "diffusion" coefficient Dt, de- 
scribes fluctuations in the rate of polymerization/depoly- 
merization. Eq. (O is identical to the diffusion dynam- 
ics of a particle in a linear potential with time-dependent 
slope —ft whose time-dependent diffusion coefficient is Dt 
(see fig. El . The novelty here is that due to the coupling 
between the MWD and monomer, both the slope and dif- 
fusivity are functionals of (j) itself. 

The reflecting boundary condition in eq. © at TV = 
represents the fact that when a living chain looses all 
of its monomers and becomes a free initiator molecule 
(i. e. reaches length A'^ = 0) it cannot depolymerize fur- 
ther and must grow again. Now we restate eq. lO as 



dt 



-vt 



dN 



D.. 



92 



0t(O)M(O) = AM , (5) 



where = d<f)t/dN. 

Now it might at first seem that despite the coupling 
of eq. (O, the monomer subsystem would effectively un- 
couple from the MWD dynamics, since one might naively 
expect that the monomer dynamics do not depend on 
the shape of the polymer MWD but only on the num- 
ber of living chains. This would imply simple first order 



kinetics for m* resulting in exponential relaxation of the 
monomer concentration, independent of the MWD dy- 
namics. The coupling however arises from the existence 
of special chains in the MWD of zero length, i. e. free ini- 
tiator molecules of which there are (^t(O), which unlike all 
other chains cannot depolymerize. Indeed, the dynamics 
obeyed by mt, or equivalently Vt, are 



dvt 



Tfast 



A0t(O) 

Tfast 



Tfast = 



/ 



(6) 

We see clearly that the only aspect of the living MWD 
which monomers see is 0t(O). Eq. (jSJ is derived by calcu- 
lating dNt/dt by multiplying eq. © by N, integrating over 
all N , using the reflecting boundary condition, and then 
using eq. The relationship between Tfast and N^o fol- 
lows from eq. (jHJ below. Here we have introduced the two 
basic independent dimensionless parameters of the sys- 
tem: 

f = TT , r = . (7) 

k+mtot mtot 

The physical meaning of / is the following: if the sys- 
tem were a pure solution of unpolymerized monomer (plus 
solvent and initiator) / would be the ratio of backward 
to forward polymerization velocities. The value of / is 
temperature dependent, being unity at the polymeriza- 
tion temperature and smaller or larger than unity in the 
polymer and non-polymer phase, respectively. The param- 
eter r, namely the ratio of living chain to total monomer 
concentration, is independent of temperature and is al- 
ways much smaller than unity. Its smallness is related to 
the mean degree of polymerization being much larger than 
unity in equilibrium (see eq. © below). 

Now setting the time derivative in eq. JSJ to zero and 
using eq. (0) it is shown in Appendix A that in equilibrium 
the Flory-Schultz distribution of eq. is recovered and 
one has to leading order in r: 



l-f 



mtot 



f 



V 



D^ 



(8) 
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Fig. 5. Schematic of tiie MWD dynamics as analyzed by 
Miyake and Stockmayer |44f and their theoretical predictions 
regarding the timescales involved, (a) Initially {t — 0) all ini- 
tiators are unpolymerized. (b) For t <^ rfast chains grow coher- 
ently with a narrow MWD until its peak reaches the equilib- 
rium mean length Naa at times of order Tfast • (c) During a third 
stage lasting up to Taiow, which was beyond their theoretical 
analysis, the MWD spreads to the equilibrium MWD. 



where r ^ 1 and subscript oo denotes the t ~* oo equi- 
librium value for the corresponding variable. We assumed 
(1 — f)lr^^'^ ^ 1, i.e. that the temperature is not ex- 
tremely close to the polymerization temperature. 

An important feature in eq. JHl is that in equilibrium 
the velocity, i. e. the slope of the "potential" term in eq. 
JSJ, settles down to a very small negative value. Were there 
no diffusion, living chains subject to a negative velocity 
field would shrink to zero length. However due to the small 
magnitude of the velocity, diffusion is strong enough to 
broaden the MWD up to distances of order Nqo-, as illus- 
trated in fig. 121 We show in the following section that any 
apparently small external perturbation has an enormous 
effect on the velocity which, depending on the direction 
of the perturbation, may become either very negative or 
very positive. 



3 Response to T-jump; Positive Velocity 
Boost 

Perhaps the simplest way to perturb an equilibrium liv- 
ing polymer system is by a small temperature change (T- 
jump). Taking as an example a-methylstyrene, the data of 
fig. show equilibrium fraction of polymerized monomer 
as a function of temperature; evidently, changes in T lead 
to different values of the equilibrium monomer concentra- 
tion moo and therefore of as well. In this section we 
consider relaxation dynamics after an equilibrated system 
is subjected at t = to a small temperature change 8T 
such that the system will eventually reach new equilibrium 
values TOoo and TVqo. That is, we follow the dynamics of 
the transition from an old equilibrium towards a slightly 
different new equilibrium state. 

The magnitude of the perturbation is measured by the 
small parameter 

6^^, (9) 



^(N) 




Fig. 6. Living polymer dynamics are equivalent to diffusion in 
a potential with slope —vt- In equilibrium the slope reaches a 
very small value of order 1/iVoo (solid line) resulting in a broad 
MWD. Under a small T-jump the slope is strongly perturbed 
and depending on the direction of the jump becomes either 
positive (dashed line) or negative (dotted line). 



The value of e is related to the magnitude of the T- 
jump. For example in the case of a-methylstyrene, assum- 
ing that the system is initially below the polymerization 
temperature Tp in the region of fig. [3 between 260° K and 
280°K, one has using eq. Q 



(5T 
'50°K 



(a-methylstyrene) 



(10) 



Thus for a-methylstyrene a temperature increase, bT > 0, 
results in a negative ^ttiq, i.e. a reduced initial monomer 
concentration relative to the destination equilibrium value. 
Now the i = relative perturbation of the MWD is 



This is of order e and smaller than unity for all N 
from the unimportant values N 3> Nqo where the 
is exponentially small) as shown in fig. IHl^a). 

Because velocity, monomer concentration and 
chain length are all linearly related (see eqs. (01, 1 
Q) changes in these quantities are simply related 



1-/ 
/ 

(11) 
(apart 
MWD 



mean 
1) and 



as 



Smt 



SNt 



Svt 



(12) 



These relations are true for all times. They allow us to 
follow the dynamics of velocity alone. Once vt is known, 
Nf, mt, Dt are determined. 

From eqs. H12|l and the initial relative changes in 
V, D are: 



5vq 



SDq 



(13) 



where we define 5mt = rrit — rrioa and similarly for other 
quantities. All equilibrium values refer to the destination 
{t — oo) equilibrium values. Thus Smo is the initial (t = 0) 
deviation from the final equilibrium. 



The important point is that since iVoo ^ 1, the relative 
change in v is much larger than e : the velocity is highly 
perturbed, as shown in fig. El Depending on the sign of e, 
the velocity may remain negative as in equilibrium, or it 
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Fig. 8. (a) Sketch of the initial and final MWDs during a T- 
jump experiment. Both distributions are exponential. Though 
very close to one another, in the course of relaxation the time- 
dependent MWD becomes very different to the final exponen- 
tial MWD. (b) General perturbation (non-exponential initial 
MWD). Most T-jump results derived here apply to this general 
case also. 



250 260 270 280 290 300 

tCk) 



Fig. 7. Fraction of polymerized monomer in equilibrium as a 
function of temperature for living a-methylstyrene in tetrahy- 
drofuran solvent, initiated by sodium napthalene. Data from 
ref. Ratio of mole fraction of total monomer, mtot, to ini- 
tiator is 0.0044. Mole fraction of total monomer is 0.1538. The 
polymerization temperature Tp is near 295°K. 



may change sign. Since the relative perturbation in D is by 
contrast small, we see that the delicate velocity-diffusion 
balance which sustained equilibrium is now destroyed. 

We will show below that relaxation to equilibrium now 
occurs in three stages. During the first, dominated by 
translational motion, the MWD is boosted far from its 
equilibrium shape. The next stage involves diffusive restora- 
tion of the region of the MWD which suffered maximum 
distortion during the first stage. The third and final episode 
entails a very slow diffusion-dominated recovery of the 
global MWD. For the remainder of this section we treat 
the case of a positive initial velocity boost, where chains 
initially grow {6T < for a-methylstyrene). Negative ve- 
locity boosts, where chains initially shrink, are considered 
in section 5. 



3.1 Coherent Chain Growth: Hole Development 

During the first stage of the relaxation process, velocity 
dominates over diffusion since v has been so strongly per- 
turbed. Chains thus grow coherently, consuming the ex- 
cess monomer in a timescale rfast by which time transla- 
tion will have halted and a highly non-linear hole will have 
developed in the MWD (see fig.EJ. Here Tfast, defined in 
eq. © , is the time for the MWD to translate distance SNq 
and reach the destination mean, i.e. Tfast ~ ^Nq/Svo (see 
eqs. CSIl and using eq. (jHJ). 

To see all of this more quantitatively, consider the ve- 
locity dynamics eq. JBJ. Initially, the velocity term on the 
rhs is much larger in magnitude than the diffusion term. 
We show in Appendix B that this remains true up until a 
time Tqs, defined below. It follows that for these times Vt 



relaxes exponentially. The same is true of mt,Nt and Dt 
which we recall are linearly coupled to Vt (see eq. ((T^ ). 
Thus 



Svt 



ev e 



(14) 

Note that for t ^ Tfast there are no free initiators (at 
much longer times these are restored, see below). Thus 
for these intermediate times all chains are identical as far 
as on and off monomer kinetics are concerned, and the net 
rate of monomer-polymer mass exchange can only vanish 
if Vt vanishes. For this reason vt decays to zero, though 
much later it will recover its small negative equilibrium 
value. The 4>t dynamics, eq. (jSJ, thus simplify to 



dt 



o(Pt 
ON 



(i«Tqs) (15) 



where the diffusion coefficient may be approximated by 
Doo since the contribution of 6Dt is negligible (see Ap- 
pendix B). This linear equation, plus the time-dependent 
boundary condition of eq. Q , has solution 



dN'G'r''{N,N')MN') (i«Tqs) (16) 



which describes a translating and simultaneously broad- 
ening MWD, as shown in fig.El Here Gf'''''{N,N') is the 
propagator of eq. H15() whose properties are calculated in 
Appendix B. 

An important quantity in what follows is (/'t(O), i. e. the 
concentration of free initiators. Using eq. p6|l we show in 
Appendix B that 



0t(O)/</)oo(O) « 

F{t*/tfl^e-*/'' {t* «i«Tfast) 



4D„ 



(17) 



where C = 4/7r^/^ and F ^ tt"^/^. Eq. I^Tj) has a clear 
physical meaning: as the MWD translates, the position 
of its peak at A'' = after time t has shifted to chain 
lengths of order vot. Were diffusion absent, this would have 
created a depletion region in the MWD at small chain 




N„ Noo No Noo No Nco 

(a) t«t* (b) f*« t « Tfo,t (c) ^ = 'Cfast 



Fig. 9. Initial response of tlie MWD after a positive boost, case of strong perturbation (t* <^ ifast)- (a) For small times, chain 
end "diffusion" smooths out the hole which would otherwise form at small due to chain growth, (b) At larger times coherent 
chain growth beats diffusion and a hole develops in the MWD. (c) Coherent growth stops when mean chain length reaches a 
value very close to Noo, displaced by an amount eNoo/d from the original value TVq. For typical cases, the parameter 6 is of 
order unity. 



lengths. However diffusion smooths out inhomogeneities 
on distances of size {Doct)^^^- Hence for times shorter 
than the crossover time t* for which vot* ~ {Doot*y^^ , 
diffusion has enough time to fill the translationally in- 
duced hole and the relative deviation from equilibrium is 
small. For times longer than t* however, the MWD trans- 
lates a distance away from the origin much larger than 
what diffusion could have homogenized, a hole develops 
in the MWD, and the concentration of zero length chains 
becomes exponentially small. 

In eq. H17|) we assumed that the perturbation is large 
enough such that t* <C Tfast • The special case of extremely 
small perturbations where this is no longer true is dis- 
cussed separately in section 5. 

3.2 Diffusive Length Relaxation: Hole Filling 

For times t 3> Tfast the fast variables monomer and mean 
chain length have relaxed very close to their equilibrium 
values. We have seen that velocity becomes exponentially 
small. Meanwhile (see below) the number of free initiators 
4>tiO) is gradually recovering. At a certain moment, there- 
fore, the 2 terms on rhs of eq. © become equal to one an- 
other and we show in Appendix B that for all later times 
the velocity time derivative is much smaller. Hence the 
velocity dynamics are now fundamentally changed. The 
new regime is one of quasi-static evolution, enslaved to 
the dynamics of 0t(O): 

Svt^ -V'S(f>t{0) (t>Tqs) . (18) 

Physically, this reflects the fact that the the only aspect of 
the MWD shape seen by monomers is the amount of spe- 
cial non-depolymerizing zero-length chains. We see that 
as the MWD slowly rearranges itself, so the fast variables 
vt , rrit and Nt variables respond quasi-statically. 

Notice that the quasi-static regime does not onset af- 
ter Tfast but rather after Tqs. The reason is that during the 
initial boost the term DtiptiO) in eq. © decays exponen- 
tially on a timescale t*; its magnitude therefore at Tfast 
is much smaller than vt which decays on a timescale Tfast 



(see eqs. H14|) and (I17|| V There is thus a cross-over period 
before the quasi-static balance between the two terms is 
established. 

Now let us examine the MWD dynamics. Since vt be- 
comes exponentially small after the MWD stops translat- 
ing, we have 

d(hf d'^chf , X , s 

^ « i^oo ^ (Tfast « t « Tfill) (19) 

with reflecting boundary conditions at the origin. Here 
once again the contribution of &Dt is negligible (see eqs. 
(I12|l . 118|l . and (|21ll ) and the corresponding term has been 
neglected in eq. (|19|) . Here 

Tfill = ^ Tslow , Tslow = . ^ ■ (20) 

The timescale TsIow is the longest relaxation time of the 
system, the time for the slow global shape characteristics 
of the MWD to relax. It equals the diffusion time for the 
MWD width Noo- Meanwhile Tfm is the diffusion time for 
the hole width, eNoo/0- It is shown in Appendix B that 
eq. H19|l has solution 

1 /2 

^ H e-^«"/* (Tfast « t « Tfiii) (21) 

where H — vr^^/^. Thus for times shorter than Tfin the 
concentration of zero length chains remains exponentially 
small while for longer times the hole fills (see fig. I10|l and 
(/)f(0) recovers its equilibrium value. 



3.3 Linearized Dynamics at Long Times 

For times beyond Tfin the hole-filling is complete and the 
MWD's non- linear feature has disappeared. Thus, finally, 
a truly linear regime onsets: relative deviations of all vari- 
ables from equilibrium are less than unity and perturba- 
tion theory can now be applied. This is done in Appendix 
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C where we find 

^ « -SMO) « 

V 

J0oo(O) (TfinA)l/2 

Me0oo(O)(t/T,iow)- 



3/2g-t/(16T,io 



lo 

{t > "^slow) 



(22) 



where J, M are positive constants of order unity. Thus 
(50t(O), and all of the fast variables which remain quasi- 
statically enslaved to 0t(O), relax exponentially on a time- 
scale Tsiow with a power law prefactor. 



4 The Second Moment of the MWD Relaxes 
More Slowly than the First 

A picture has emerged of fast mass- related variables which 
relax in Ttast and slow MWD shape properties which relax 
in the much longer time Tgiow The MWD's first moment 
defines the total mass in the polymer system: this is a fast 
variable and was considered in section 3. The simplest slow 
shape property is the dispersion At, 



At EE Ni 



dNN^MN) , (23) 



which is closely related to the second moment. In this 
section we consider the very different and slow relaxation 
of this quantity. 

Substituting eq. Hll|) in eq. (|23|l one finds that fol- 
lowing a T-jump the initial relative perturbation in the 
dispersion, SAq, is 



^^0 _ g £ 
A^ 9 ' 



(24) 



where the equilibrium value, Z\oo, is found using eq. (|T)). 
Now the dynamics of the dispersion can be derived by 
evaluating the time derivatives of N^, Nt in eq. (|^ . These 
can in turn be calculated by multiplying eq. ^ by N'^ 
and N, respectively, integrating over all N and using the 
boundary condition of eq. JSJ- One has 



1 d6At 
dt 



2r. 



low 



Do 



Nt 0t(O) 

A^oo0oo(O) 



1 ^0t(O) 



2tsi, 



ow V^ctt 



(0) 



(25) 



Here we replaced Dt and Nt with their t = oo values. This 
is correct to leading order since the relative perturbations 
in these quantities are much smaller than that of '?!>t(0) 
which undergoes much larger changes (see section 3). 

Now integrating eq. 1(23) using eqs. (ETj), ijllj and if^ 
one finds that up to the time Tf^\\ the dispersion has changed 
very little relative to its value just after the perturbation 
at t = 0: 



At-Ap 
A^, 



r(C/3)(t/r3iow)(tA*) 

\t/(2Tslow) 



1/2 



{t < t*) 

{t* -^t-t: Tfill) 



(26) 



^(N) 



m- ( 




Fig. 10. Schematic of MWD during hole-filling episode. Dot- 
ted (dashed) lines show MWD profile for times rfast ^ i ^ Tmi 

(Tfill < t < Tslow). 



This is in stark contrast to the first moment which had 
already relaxed by the much shorter time Tfast- The dis- 
persion's relaxation process does not properly get going 
until times of order TsIow: 

SA ^ ( SAo/A^ + Je6'-i(t/T,iow)^/^ (Tfiii < t < t,iow) 
A^ ~ \ae (VT-siow)-'/'e-*/(if^"='-) [t » t^iow) 

(27) 

where a = 8A//39. The relative perturbation in the dis- 
persion remains of order e for t <C Tgiow and subsequently 
relaxes exponentially. 



5 Response to T-jump; Negative Velocity 
Boost 

In this section we study the case opposite to section 3, 
when e < and the T-jump induces a higher equilibrium 
monomer concentration and a smaller N^c (see eq. ©)• 
Such a perturbation in the example of a-methylstyrene 
would be produced by an increase in temperature as one 
may see in fig. [7| In this case the sign of the initial per- 
turbed velocity (eq. (|13|l ) is negative (see fig. (HI). Thus 
initially the MWD is not boosted towards larger molecu- 
lar weights as in section 3, but is instead boosted in the 
opposite direction, towards smaller N. Therefore, instead 
of depletion, an excess of small length chains and free ini- 
tiators is produced, creating instead of a hole, a sharp 
peak in the MWD near = as shown in fig. ^2 

Despite the differences between the two cases, many 
of the results of section 3 remain unchanged. In fact all 
equations of section 3 up to eq. (|16|l remain unchanged, 
the only difference being that e is negative. Thus eq. (|16|l 
now describes a MWD coherently shrinking and simul- 
taneously broadening with diffusion coefficient Doo- But 
since free initiators cannot depolymerize, excess polymer 
must build up near N = 0. Indeed, starting from eq. 
it is shown in Appendix D that 



>o(0) 



F' [t/t*) {t* <^t^ Tfast) 



(28) 
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Noo Nn 




No, No 




No, Nn 



(a) t «t' 



(b) r«i«Tj 



fast 



fast 



Fig. 11. Initial response of MWD after a negative boost (strong perturbation, t* <C Tfast)- (a) Chains shrink coherently and 
MWD moves from the initial mean iVo towards iVoo- For small times, diffusion smooths out excess peak at small A'^. (b) Beyond 
t* diffusion is too slow to counteract coherent chain shrinkage and short chains accumulate at the origin, (c) Coherent shrinkage 
halts when mean chain length approaches Noo- 



where C" = 4/71^/^ and F' = 4. Eq. has a similar in- 
terpretation to eq. H17|l . For t <^t* diffusion is fast enough 
to smooth out the excess polymer accumulated by the neg- 
ative velocity at the origin and thus 0t(O) remains close 
to its initial value (see fig. ^2 (a))- For t ^ t* however the 
flow towards the origin is faster than what diffusion can 
smear out and a peak forms whose height increases with 
time, as shown in fig. ITIT b). (Once again it is assumed 
that the perturbation is not so small that t* <C Tfast; the 
opposite case is examined in section 5). 

After the completion of the first stage at times of order 
Tfast) similarly to section 3, the fast m^, Nt variables have 
essentially relaxed and start to respond quasi-statically 
to the slow diffusion-driven shape changes of the MWD. 
In fact we show in Appendix D that eqs. H14() and H18|l 
still apply, but now with a different crossover timescale, 
Tqs = Tfast ln ^oo ■ This difference arises because unlike the 
positive boost case (/'t(O) does not become exponentially 
small and so needs less time to catch up with the vt term 
in the velocity dynamics, eq. JBJ. 

Similarly to the positive boost, eq. (|19|l still applies 
for t < Tfiii. By this time the accumulation at the origin 
is able to diffuse to distances of the same size as the re- 
gion from which it was transferred from during the first 
stage. Hence (f>t{0) becomes of order its initial value, and 
thereafter perturbation theory is valid. In Appendix D we 
show that 



6t(O)/0oo(O) « 

1 + J'(Tfiu/t)i/2 
H-M'e(t/Tsiow) 



-3/2 -t/(16r,io„) 



(Tfast < t <C Tfiii) 
(Tfill < t < Tslow) 
{t > Tslow) 

(29) 



where K',J', and M' are positive constants of order 
unity. Note there is apparently no smooth cross-over at 
t = Tfast between the forms of eqs. (|28|l and H29|l. This 



is because there is in fact an extra episode at times of 
order Tfast during which the very thin peak spreads from 
width {Doot*Y^'^ to a width (ZJcxjTfast)^/^- Note also that 
the time behavior for t 3> Tfm has been derived using the 
linearized dynamics of Appendix C. 

Finally, to obtain the dispersion dynamics we use eqs. 
and (|2ni) in eq. to obtain 



SAt/A^ - SAo/Aoo ~ 

-(C73)(t/Tsiow)(t/t*)i/' 
-F74(t/Tsiow)(Vi*) 

-SAo/A^ - (8Af 739) e (ts. 



{t < t*) 

(i* < i < Tfast) 
(Tfast < t < Tfill) 
(Tfill < t < Tslow) 

3w/t)^/2e-*/(i6-.iow) 

(t » Tslow) 



(30) 



which is very similar to eqs. (|26|) and (|27|l for the positive 
boost case. Again the dispersion relaxes on a timescale 
Tslow , much longer then the relaxation time of Nt . 



6 Very Small T-jumps: e less than its critical 
value ec 

Hitherto we implicitly assumed sufficiently strong velocity 
boosts that the timescale t* after which translational mo- 
tion dominates diffusion and a hole or peak begins to form 
is shorter than the time for full development of the hole 
or peak, namely Tfast • Now the ratio of these timescales is 



t* 



Tfast 



9I/2 



1/2 



(31) 



which defines a critical value of the perturbation param- 
eter, Ec- Clearly, the system's response must be different 
for e values below this rather small threshold. We consider 
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f(0) 



\\ 80(0) - (^) ' 




Fig. 12. Schematic of MWD development during rfast <C t >C 
Tsiow when the small chain peak spreads diffusively. Dotted 
(dashed) lines correspond to times rfast <^t<^ rmi {Tmi -^t <^ 

rslow ) • 

such very small perturbations in this section. Although a 
full hole or peak does not have time to develop, we will 
find that the response is still strong and non-linear. 

Since such weak perturbations do not involve the cre- 
ation of a complete hole or a peak in the MWD, one ex- 
pects that the timescales t* and Tm associated with the 
their creation and destruction lose their physical mean- 
ing as cross-over timescales. This is shown in Appendix 
E where, using similar arguments to the ones of sections 
3 and 5, it is proved that weak perturbations are identi- 
cal to stronger perturbations except that the regime up 
to t* is truncated at rfast and the following regime (which 
for stronger perturbations lasted till rfin) deleted. In Ap- 
pendix E it is shown that for both positive and negative 
boosts, one has for the fast variables 



Svt 



e-tM„3t (t< rfast) 
-'5<?!>t(0) (rfast <i<Ts 

low J 



(32) 



Note that that the timescale rqs does not appear since its 
magnitude lies between those of t* and Tfm. 

The evolution of the MWD (see figs.[l21lllll is similarly 
shown in Appendix E to lead to the following solution for 
0f (0) in the case of a positive initial velocity: 



l-C{t/t*y/^ (t«Tfast) 
1 /2 

1 - Q (Tfill/i) (rfast < i < Tslow) 



{vo > 0) 



while for a negative initial velocity one has instead 

hiO) f l + C'(t/r)l/2 (t« rfast) 



(33) 



1 + Q'(Tfin/t)i/2 (rfast < t < rsiow) 



{vo < 0) 

where Q, Q' are positive constants. Finally, it is straight- 
forward to show that the relaxation for for t 3> rsiow is 
exponential as for the stronger perturbations, eq. H22|l . 



The most important result of this section is that even 
though the e < ec case is "weak" as compared to sections 3 
and 4, the system's response remains large and nonlinear. 
One sees from eqs. (T^ and (p^ that at rfast when the 
deviation of 4>t{0) from equilibrium is largest, the relative 
deviation of 0t(O) from equilibrium is much larger than e. 
Therefore even in the weak case a hole or peak does form 
in the MWD but of a smaller amplitude as compared to 
the stronger perturbation case (see figs. lT^IT^ . 

Now the dispersion dynamics may be similarly evalu- 
ated using eqs. (|33|) and (|34|l in eq. (|25|l . One finds that 
At relaxes once again on a timescale rgiow ■ For the positive 
initial velocity case one finds 



SAt/Aoo - SAo/Aoo « 

r (C/3)(t/rsiow)(tA*)i/' 
I Oe0-i(</rsiow)i/2 
I -5zio/Z\oo + (8Af/39)e 



[t < rfast) 

(rfast <^t^ rsiow) 

rsiow/t)^/2e-*/(i6-=- 



(35) 



(t 3> rsiow) 
while for a negative initial boost 



\l/2 



SAt/A^ - SAq/A 

-(C'73)(t/rsiow)(t/i 
-Q'e0-Ht/rsiow)i/' 
-SAo/A^ - (8AfV39) e (rsiow 



(t < rfast) 

(rfiii < t < rsiow) 



/t) 



5/2, 



{t > rsio 



) 



7 General Perturbations 

We have so far considered a specific type of perturbation, 
namely T-jumps. In this section we consider general per- 
turbations. We will see that almost all of the phenomenol- 
ogy remains unchanged. 

Let us start with another simple type of perturba- 
tion: addition of a very small amount of monomer to an 
equilibrated system. (This would require special care to 
avoid introducing impurities which in current experimen- 
tal methods are destroyed by the initiators before poly- 
merization |2|.) Changing the total amount of monomer 
perturbs both independent system parameters, r and / 
(see eq. JTl), while a T-jump perturbed only /. Despite 
this, the results of the previous sections remain identical. 
The reason is that since in both cases the initial MWD is 
an exponential distribution whose mean is different from 
the new target equilibrium MWD, the only variable pa- 
rameterizing the perturbation is SNq, or equivalently SrriQ. 
Thus one sees from eq. Q that monomer addition is sim- 
ply equivalent to a positive velocity boost induced by a 
T-jump. 

Consider now the most general perturbation, i. e. one 
generating an initial MWD of arbitrary shape as shown in 
fig.|Sfb). (This should be compared to the T-jump and m- 
jump which produce initial MWDs of exponential shape 
with a mean slightly displaced from equilibrium.) Such a 
general perturbation is small provided 



\SMN)\ 

0oo(^) 



(37) 
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N„ N„ 



(a) t«Xf., 



N„ N„ 
(b) t = X,i 



Fig. 13. Weak perturbations, e < e^. For such small perturba- 
tions t* ^ Tfast- (a) The coherent component of chain growth 
shifts the MWD mean towards Noo- Diffusion is dominant, 
smoothing depletion at small A^. The net effect is a nonlin- 
ear deviation from equilibrium, (b) By rfast coherent growth 
halts, when mean chain length approaches N^o- The behavior 
for longer times, t S> Tfast, is essentially the same as that shown 
by the dashed lines in fig. ^| 



for all N . As for the T-jump case, we have now 



5mQ 



5No^ 



N, 



oo Jo 



dNNSMN) (38) 



Since eq. 1)38(1 is of the same form as eq. all the argu- 
ments of section 3 leading to the large relative perturba- 
tion in V (eq. ((13|l ) remain unchanged. Moreover, in our 
analysis of sections 3, 4, and 5, the exact shape of the per- 
turbed MWD did not matter in any of the leading order 
results we obtained. The reason is that the effect of the 
velocity field on the MWD is so drastic that to first order 
the response of any MWD close enough to equilibrium is 
the same, in the sense that the MWD simply translates, 
creating a hole (or peak). Thus the analysis of sections 3, 
4 and 5 directly apply to the general case. 

A special exception is when the first moment of the 
perturbation S(f)t{0) is arranged to have a value much 
smaller than eNoo- In this case the relative amount of 
monomer-polymer mass transfer is much less than e and 
this in turn reduces the velocity perturbation. For exam- 
ple, if the first moment were chosen to be so small that 
^mo/moo ~ l/A"oo, then from eq. H13|l one sees that the 
relative perturbation in vt is of order e, much less than the 
order eiVoo produced by a T-jump. In such cases, linear re- 
sponse applies for all times. However, if the first moment 
of S(f>o is much greater than unity the response remains 
nonlinear. 



8 Discussion 

The peculiarity of living polymers is that different mo- 
ments of the MWD relax to equilibrium on different time- 
scales. In the present work we have shown that the first 
moment, i.e. the mean length, Nt, relaxes on a timescale 
Tfast- However the shape of the MWD measured by the 
second moment, or equivalently the dispersion. At, relaxes 
after TsIow Because the two timescales depend on different 



,1 5(|)(0)~(f/f*j" 




N„ No 



(a) t«Xf, 



N„ No 
(b) t = Xf,„ 



Fig. 14. As in fig. ITsl but now for a negative initial velocity 
boost; a peak is formed instead of depletion. Evolution for 
t 2> Tfast as dashed lines of fig. fT^ 



powers of Noc their ratio may be extremely large: 



Tfast 



/ 



V- (l-f) ' 



Tslo 



iV2 

4v- 



(39) 



where l/v~ is the natural time scale and / « "^oo/"^tot is 
approximately the fraction of non-polymerized monomer. 
Our only assumption has been that / is not extremely 
close to unity, i. e. the system is not very close to the 
polymerization temperature. 

In addition to the MWD, the other important observ- 
able is the free monomer concentration, mt- We emphasize 
that since the number of chains is fixed by the number of 
initiators, thus mt is linearly coupled to Nt by mass con- 
servation for all times. Hence the relaxation of these two 
quantities is identical and mt relaxes after Tfast • Note that 
for certain other living polymers such as end-polymerizing 
wormlike surfactant micelles, this is not true |29I3(JI5!3] . In 
these systems there is no separate initiator species and a 
pair of monomers may spontaneously unite to form a living 
chain. The dynamics are thus very different and mt and Nt 
may relax on different timescales. Marques et al. stud- 
ied this class of living polymer; after linearizing dynamical 
equations they found that after a small T-jump mt relaxes 
much more rapidly than Nt which subsequently relaxes in 
time ~ N^. In numerical simulations of this same class 
of system, Milchev et al. PHI report a, 1/t decay of mean 
chain length after a time ^ N^. In the related system 
of spherical micelle aggregation |i56>.57. .58 , 59 , the num- 
ber of micelles can change and again the kinetics are very 
different to the present case. 

A summary of our findings is shown in fig. 1151 for the 
case emphasized here where a small T-jump induces a pos- 
itive velocity boost to chain growth. The figure shows the 
very different relaxation times of first and second MWD 
moments. It also shows the interesting behavior of the 
number of free initiators, ^t(O), which suffers enormous 
depletion at intermediate times despite ultimately recover- 
ing to a level very close to its initial value. We showed that 
whilst virtually all of the relaxation of the fast variables 
Nt,mt occurs on a timescale Tfast, they are actually not 
completely relaxed by this time. The final very late stages 
of their relaxation during which their values are fine-tuned 
(to relative order 1/iVoo) to the final equilibrium values. 
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Fig. 15. Time dependent recovery following a small perturba- 
tion, positive initial boost, e > ec. Top: SNt, Smt (solid line) 
become very small after Tfast while the MWD dispersion 5 At 
(dashed line) relaxes after TsIow Bottom: In the process of equi- 
libration free initiator concentration undergoes large changes 
even though its initial and equilibrium values are very close to 
one another. 



occurs on a timescale Tgiow • During this episode their val- 
ues evolve quasi-statically, enslaved to (j>t{0) according to 
6Nt9/Noo ~ —Smt/moo ~ <5(/>t(0). The novel behavior of 
^t(O), as well as its central role in late relaxation of the 
fast variables suggests this as an interesting quantity to 
measure experimentally. This might be achieved spectro- 
scopically and would have the advantage of being a rela- 
tively uninvasive probe. We remark that for bifunctional 
initiators, e.g. those used in studies by the Greer group, 
the meaning of (f>t{Q) is the concentration of "half-chains," 
i. e. those having at least one chain end which is a bare 
initiator molecule. 

We have demonstrated that the type of dynamical re- 
sponse depends on the magnitude and sign of the pertur- 
bation. For example, if the temperature change is reversed 
in sign, then the velocity boost also changes sign. Thus for 
a-methylstyrene a negative T-jump produces the positive 
velocity boost phenomenology of fig. ^| whereas a small 
temperature increase causes a negative boost. In this case, 
instead of short chains being annihilated, their number 
increases dramatically during the relaxation process. As 
far as perturbation magnitude is concerned, if the relative 

~ 1/2 

value e is less than a critical value ^ 1/Noo then the 
response is milder but remains non-linear. For a positive 
boost, the hole in the MWD is no longer complete, but 
instead has a relative depth less than unity but still much 
bigger than e. 

How do our predictions compare with experiment? In 
ref. small T-jumps were imposed on a-methylstyrene 
systems in the semi-dilute regime at concentrations high 
enough to be nearly entangled. Relaxation was monitored 
by measuring viscosity 77 as a function of time. Now gener- 
ally we expect viscosity to scale as a polymer-concentration- 
dependent power of mean chain length, rj ^ c^-ZV''', where 



the value of 7 is predicted by standard theories of polymer 
physics |6UI47j . For a general MWD, note the coefficient 
which depends on the shape of the MWD, i. e. it depends 
on the full function [^(iV)]. (This is because 77 depends not 
only on the first moment, but in general has a complex de- 
pendence on chain length distribution). Thus we predict 
an initial fast relaxation of t] lasting time Tfast — N / {9v~) 
(corresponding to the relaxation of N) followed by a much 
slower relaxation in Tgiow (this is the relaxation of the pref- 
actor C0). The fast relaxation time is independent of T- 
jump magnitude (i.e. independent of e). Now from eqs. 
(|SJ| and O we can rewrite Tfast = ^/{rk^mtot)- Using 
r — 4.7 X 10^^, mtot = Q.29gm/cm^ and a measured 
value 12 fc+ w 0.2M~^sec~^ one estimates Tfast ~ 1000 
sees. (Note that this work involved bifunctional initia- 
tors. In ref. |54| we incorrectly estimated Tfast to be twice 
the value calculated here based on a wrong estimate of 
the number of living chains which was taken to be equal 
to the number of bifunctional initiators instead of twice 
this value.) Similarly, with j2| / ~ 0.5, v~ = O.lsec"^ 
one finds TsIow ~ 1 month. Hence for these experiments 
(timescales much less than months) ij should relax in ap- 
proximately 1000 sees. In the experiments, jumps ST « 1 
"^C were imposed at various temperatures in the range 
283°K :$ T < 290°K. Now the temperature dependence 
of Tfast ~ 1/A:+ is determined by A:+ which is reported jT2] 
to vary by « 10% over this temperature range. Thus we 
expect only a slight variation in the viscosity relaxation 
time at the different temperatures studied (despite the 
fact that Noo changes significantly) . These predictions are 
very close to the experimental findings, where relaxation 
times at all temperatures were found to be about 2000 sec- 
onds, with very little variation from one temperature to 
another. This is in sharp contrast to the relative changes 
in r] itself, which varied by almost an order of magnitude. 

We have emphasized small perturbations in this work. 
In fact, the response to a large perturbation is obtained, 
qualitatively speaking, by setting e = 1 in our results. This 
case was addressed by Miyake and Stockmayer 0^ who 
considered an initial condition in which all initiators are 
free (i. e. Nq = 0). They were the first to identify the two 
timescales of the system, Tfast and TsIow, in the limit of 
very small depolymerization rates (/ <C 1). Their analyti- 
cal results applied to the t < Tfast stage where they found 
an initially uniformly translating and broadening MWD 
peak, reaching iVoo after Tfast as shown in fig. [5| Seen from 
the viewpoint of a general perturbation of magnitude e, in 
this case Tan and TsIow coincide because e = 1. That is, the 
"hole" (in this case the entire region between and iVoo) 
fills up on the same timescale in which the whole MWD 
relaxes. The linearized dynamics, which for a small pertur- 
bation appfied for t > tru due to 6(j)t{0)/4>oc,{0) becoming 
small (see section 3.3) cannot be employed. Describing dy- 
namics at times of order TsIow in this case is thus a difficult 
nonlinear problem j5 1115211^ . 

In recent experiments, Greer's group ^1^] studied 
such large perturbations using a-methylstyrene. After a 
rapid quench below the polymerization temperature, mono- 
mer concentration and MWD were analyzed after a time 
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delay. The time dependent relaxation was probed by re- 
peating the procedure for different time delays and sam- 
ples. They found different relaxation times for the mono- 
mer concentration and the MWD width, as predicted by 
theory However, Nt and rrit did not have the same 
relaxation, as mass conservation would seem to dictate. 
This may be due either to side-reactions during polymer- 
ization or living chain ionic association effects (6f..62M63| 
which have been studied theoretically Dynamic light 
scattering measurements (i 1 suggest that prior to poly- 
merization the initiators self-assemble into long polymeric 
structures and it has been suggested [12] that initiators 
may not all be equally available for polymerization follow- 
ing the T-quench. However, lifetimes of aggregate struc- 
tures would need to be extremely long for these effects to 
interfere with polymerization dynamics since as we have 
seen the MWD relaxation times are very large. 

In conclusion, our hope is that this work will motivate 
further experimental study of the dynamical sensitivity of 
living polymers to small perturbations. The drastic effect 
on the MWD in the small chain region, in particular the 
depletion or excess of free initiators, is a natural focus for 
experimental measurement. One can think of other time- 
dependent perturbations such as small amplitude thermal 
cycling which would probe interesting aspects of their ul- 
trasensitive dynamics. Finally, the systems we have ana- 
lyzed here are model starting points for the considerably 
more complex biological living polymers, actin filaments 
and microtubules, which are intimately involved in the lo- 
comotion and structural integrity of living cells |39l4l)l4l] . 



This work was supported by the Petroleum Research Fund, 
grant no. PRF-33944-AC7. 



A Derivation of Equilibrium Values moo. ^ooi 

The equilibrium MWD can be calculated by setting the 
time derivative in eq. jSJ to zero. Its solution is the Flory 
distribution, eq. where 



B Self-Consistency of Results of Subsections 
3.1 and 3.2 

B.l Coherent Chain Growth 

The solution for the Vt dynamics, eq. (|14|l . follows by solv- 
ing eq. ijnji after neglecting the (f>t{0) term on the rhs. 
For t <^ Tqs, this is a self-consistent solution since then 
the (/'t(O) term is negligible with respect to the remaining 
terms as can be seen using the solutions for vt,4>t{0) of 
eqs. inil, (UZIl, and eq. 

Eq. lfT31) is derived by substituting Vt from eq. ffT^ 
in eq. (O, replacing Dt = Doo + 5Dt, and dropping the 
5Dt term. The fact that this term can be neglected may 
be seen as follows. Even if 5Dt did not decrease but re- 
mained of order 5Dq throughout this regime this would 
lead to replacing D^o by I?oo(l + e) in all results obtained 
in section 3 (see eq. CSJ)- This upper bound on the ef- 
fect of 5Dt can easily be seen to lead to higher order in 
e terms (see eg. eq. (jlTIl ') and the 6Dt term can therefore 
be neglected. 

Now let us consider the propagator of the dynamics 
of eq. (|f 5|l . which appears in eq. H16|l . For short times, 
t ^ Tfast, one may replace Vt in eq. (|15|l by vq. Eq. (|15|l 
then has constant coefficients and its propagator including 
the boundary condition may be calculated exactly (HH] : 

I 



^^-(N-N' -votf /(AD^t) _^ ^-N'vo/D^^~(N+N'-votf/{4:D^t) 

Here erfc(a;) = 1 — erf (a;), where erf is the error function. 
The solution for 0t(O) is derived by substituting eq. 
in eq. lfTC|l and setting iV = 0. One has 



■ erfc(a;) + 2x^eYic{x) ~ 2TT-^''^xe-^^ 



1/2 



(i<Tfast) (43) 



Solving the system of eqs. 
has 



1 + / + 



(HOI) and 



© for 



(40) 



, iVoo one 




(l-/)2 + r(I + 3/) + 



1/2 > 



"^tot 



Here we used the fact that since the values of N' con- 
tributing to the integration are much smaller than iVoo, 
one may approximate 4>o{N') by (/)oo(0) in the integrand 
of eq. (|f 6(1 (with relative errors of order 1/iVoo). Eq. 1(17(1 of 
the main text follows from eq. 1(43(1 by taking the two cor- 
responding limits. Note that eq. ((43(1 approximates (/'t(O) 
for all t <C Tfast, including times of order t* . 



(41) 



For r ^ 1 and assuming (1 — f )/r^^^ ^ 1, by expanding 
eq. ((41(1 in powers of r and keeping only the first term of 
the expansion one recovers eq. ^ . The values of Woo , D^o 
in eq. ^ are obtained by substitution of moo , -^oo in eq. 



B.2 Hole Filling 

Eq. ((18(1 is a self-consistent solution of eq. ^ since for 
t 3> Tqs, assuming the validity of eqs. ((18(1 and ((17(1 one 
may verify that eq. jnj is satisfied with the term on the 
Ihs being much smaller in magnitude than both terms on 
the rhs which balance one another. 
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Consider now the validity of eq. ((11^ for Tfast ^ t <C 
Tfiii. Substituting eq. (|21|l into eq. H18(l one has explicit 
solutions for rrit, Nt,Vt, Dt- One thus sees that in the time 
regime of consideration the velocity term in eq. © is ex- 
ponentially small and therefore can be deleted in eq. H19|l 
with very small error. 

The MWD dynamics are therefore described by eq. 
H19|) whose propagator satisfying reflecting boundary con- 
ditions is 



Gf"(7V,iV') = 



1 



|g-(W-W')V(4Doot) _,_ g-(JV+JV')V(4D„„t)| , (44) 

The solution of eq. 119|l is thus 



diV>*„_(7V')G^l_(iV,iV') 



dN' MN' - eN^/e)Gf'''{N,N') 

(Tfast < t < Tfiii) . (45) 

where Across Satisfies Tfast < Across < Tqs, i. 6. it belongs 
to the time regime described simultaneously by eqs. H15|l 
and (|19|l . In going from the first integral in eq. l|l5l) to 
the second we used the fact that at tcross the MWD has 
approximately the shape of the initial MWD, boosted in 
the positive direction by eNao/0. The effect of diffusion 
during the initial boost of the MWD may be shown to be 
small. Notice also that evolution in the last expression in 
eq. (|45|l starts at t = 0; for t ^ icross we may to leading 
order replace t — Across in C*'^ in eq. by t. 

Substituting eq. (|44|l in eq. (|45() and setting = one 
obtains 



«(0) 



rfc 



Cf'f 



(Tfast < t < Tfill) (46) 



Considering times much shorter than the crossover time 
Tfill, eq. is recovered. 



C Linearized Dynamics, t ^ Tfiu 

Ultrasensitivity is such that the dynamics cannot be lin- 
earized until the later stages, t 3> Tfin. For these times we 
may linearize eq. lO by dropping terms proportional to 
products of (5(/>t, 6vt, SDf. One has 



dt 



ON 



6vt 



(47) 



with reflecting boundary conditions at the origin. In eq. 
(|47() the 5Dt term was neglected since its magnitude is 
smaller by l/N^o than the corresponding Svt term. The 
(5-function in the source term in eq. H47|) arises from lin- 
earizing the boundary condition in eq. ij^l. Notice that 



by definition, (5(/'t(0) is normalized to zero and that the 
source term in eq. (|47ll preserves this normalization since 
its integral over all N is zero. 

Performing the corresponding linearization in eq. © 
one has 



-rSvt 
dt 



1 



1 



Svt DooSMO) 



Tfast 



Tfast 



(48) 



Now since eq. (|47|l is of the same form as eq. (|15|) . 

its propagator, G\^^°^'^{N, N'), is given by eq. ig^l with 
Vq replaced by v^o- Thus if (5(/>tj^ (N) is known at time i^,, 
then the solution of eq. (|47|l for subsequent times is 



SMN) « / dN'G'rtT{N,N')5c^tAN') + Rt{N) , 
Jo 

(49) 

where 



Rt{N)= dN' dt'G';^1f'{N,N')fit'{N') 



' dt' [0^(jV)-Gl'— (A^,0)] • (50) 



Here the last expression for Rt (N) is obtained by perform- 
ing the N' integration using the expression for /if from eq. 
(I47|l and the fact that (poo{N) remains unchanged when 
evolved with G""'=^''. In eq. ^ we also used 



Svt ~ -D^dMO) 



(51) 



which is valid throughout the linearized time regime as 
one may derive using eq. H48|) . The solution of eqs. (|47|l 
and (|48|l is thus found by setting = in eq. 1)49(1 , solving 
for SiptiO), and then using the calculated expression in eqs. 
(gni and (EU to find <5<^t(A^) and Svt- Let us now perform 
this analysis for two time regimes. 

(i) Tfill ^ i ^ Tgiow In this regime, we will see the 
source term can in effect be ignored. The top expression 
on the rhs of eq. (|22() is a self-consistent solution of eq. (|49|l 
which is proved as follows. As one may see from eqs. (|21|l 
and (|22|) there exists a constant f3 of order unity such that 
for t > /Sthii the relative perturbation in iptiN) is much 
smaller than unity and eq. l|T7)l applies. (More precisely, 
for any desired relative smallness there exists a different 
constant /3.) Thus setting ti, = (3Tf^\\ and A^ = in eq. 
and using the expression for 5<t)t (0) from eq. 1(22(1 one finds 
that with appropriate choice of J, eq. 1(49(1 is satisfied with 
Rt (0) being much smaller than the remaining two terms in 
eq. ((49(1 . The integral term in eq. 1(49(1 is calculated using 
the fact that as seen from eqs. ((21II . (I22II . and fig. 
5(t)tdN) is of order (5(/>tjO) > e<Poo{N) for A^ < eN^/O 
and of order e(j)oo{N) for larger A^. We also used 



Gn'"(0,iV') 



l/[D^{t-t')f' {N' ^ D^{t - t')) 
{N' :^ Doo{t - t')) 

(t,t'<Tslow) (52) 



which may be derived from eq. ((42(1 . 
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(ii) t » Tsiow Using eq. (|^ one finds 



linear 



G 



where 



(0,iV') «5oo + (5t-5oo)A(7V') (i»r,iow), (53) 



St = G'r^'iO, 0) , Soo= cl^ociO) = l/^oo , 



X{N') 



N' 
2nZ 



(54) 



Now setting II — TsIow in eq. I^?^ and using eqs. 
and one has 



bt{Q) ^{St-S^) / dN' \{N')5^r^,^^{N') + 

JQ 



1 



45', 



oo'^slow JTsI, 



dt' {St-f -S^)6cPt,{0) (55) 



where we used the fact that S(I)t^-^^^{N') is normahzed to 
zero. 

Now the self-consistency of the t > Tgiow expression on 
the rhs of eq. (|22(l is proved by substituting eq. (|22(l in eq. 
H55|) . Since Scpr^ia^iN') is of order e/Noc, the magnitude 
of the integral of S(j)r^^^^ {^') in eq. 1)55(1 is of order e. The 
last integral term in eq. ((55|l is evaluated using eq. H22|l 
and 

St >S'oo ~ 

165oo7r-i/2 (T,iow/i)'/' e-*/(i6-=.ow) » ^^j^^) 

(56) 

which may be proved using eqs. (|^ and One finds 
that for t 3> Tgiow all terms in eq. (|55|) have the same 
time dependence and are of the same order of magnitude. 
Thus eq. 1(55(1 may be satisfied by appropriate choice of 
the numerical coefficient M in eq. ((22(1 . 



D Self-Consistency of the Results of Section 
5 

D.l Coherent Chain Shrinking 

The validity of eq. ((14(1 is verified using the same argu- 
ments as in the first paragraph of Appendix B but now 
using the expression for (j>t (0) from eq. I(28() . 

Setting iV = in eq. one obtains eq. (|^ with x 
now replaced by — x. Considering times much greater and 
much less than t* leads to eq. ((28(1 . 



D.2 Peak Decay 

The validity of eq. I(18() is proved using exactly the same 
arguments as the ones in the first paragraph of the corre- 
sponding section of Appendix B, but now using the expres- 
sion for (ptiO) from eq. ((291) . One finds that the crossover 



time Tqs ~ Tfast lii(iVoo) is different from the respective 
crossover time found in Appendix B. 
The solution of eq. 1(19(1 is now 



MN)- / dN'MN' -eN^/0)Gf*'{N,N') + 
Jo 

Gf"" {N, 0) / dN' 00 {N') (Tfast « t « r,iow) (57) 

Jo 

Here the evolved MWD consists of two parts, correspond- 
ing to the two terms on the rhs (see fig. Illf c)). The first 
is the initial MWD, (po, shifted by eNoo/d towards smaller 
N . The second part is the excess peak accumulated by 
the negative boost at the origin and whose amount was 
the polymer initially lying between zero and eNoo/0. Here, 
as in section 3, the effects of diffusive broadening during 
t < Tfast have been neglected, and Tfast in eq. 1(57(1 has been 
replaced by zero. 

Now substituting eq. 1(44(1 in eq. ((57(1 and setting N — 
one obtains eq. ((29() similarly to the derivation of eq. ((21(1 
in Appendix B. 



E Self-Consistency of the Results of Section 6 

For times t <^ Tfast, all results and analysis are identi- 
cal to those for the strong perturbation (e > Cc) case in 
the regime t ^ t* . Eq. 1(32(1 can be shown to be a self- 
consistent solution using the same arguments as those of 
Appendices B and D to prove eqs. ifT^ and ifT^ . respec- 
tively. For times longer than Tfast hncarization of the dy- 
namics can be performed and the results of Appendix C 
directly apply. 
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Abstract. We study theoretically the dynamics of living polymers which can add and subtract monomer 
units at their live chain ends. The classic example is ionic living polymerization. In equilibrium, a delicate 
balance is maintained in which each initiated chain has a very small negative average growth rate ( "veloc- 
ity") just sufficient to negate the effect of growth rate fluctuations. This leads to an exponential molecular 
weight distribution (MWD) with mean N. After a small perturbation of relative amplitude e, e.g. a small 
temperature jump, this balance is destroyed: the velocity acquires a boost greatly exceeding its tiny equi- 
librium value. For e > ec ~ l/N^^^ the response has 3 stages: (1) Coherent chain growth or shrinkage, 
leaving a highly non-linear hole or peak in the MWD at small chain lengths. During this episode, lasting 
time Tfast ~ .^V, the MWD's first moment and monomer concentration m relax very close to equilibrium. 
(2) Hole-filling (or peak decay) after tah ~ e^N^. The absence or surfeit of small chains is erased. (3) 
Global MWD shape relaxation after TsIow ~ -^V^. By this time second and higher MWD moments have 
relaxed. During episodes (2) and (3) the fast variables [N , m) are enslaved to the slowly varying number 
of free initiators (chains of zero length) . Thus fast variables are quasi-statically fine-tuned to equilibrium. 
The outstanding feature of these dynamics is their ultrasensitivity: despite the perturbation's linearity, the 
response is non-linear until the late episode (3). For very small perturbations, e < Ec, response remains 
non- linear but with a less dramatic peak or hole during episode (1). Our predictions are in agreement with 
viscosity measurements on the most widely studied system, Q-methylstyrene. 

PACS. 82.35.-X Polymers: properties; reactions; polymerization. 05. 40. -a Fluctuation phenomena, random 
processes, noise, and Brownian Motion. 87.15.Rn Biomolecules: structure and physical properties; Reactions 
and kinetics; polymerization. 
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Living polymers are intriguing examples of soft matter spond dynamically via polymerization and depolymeriza- 

and have major roles in technology and biology. The clas- tion reactions allowing a new thermodynamic equilibrium 

sic system, of major importance for the synthesis of high- to be attained. These systems are "living" in the sense 

grade polymer materials and the subject of a long his- that a change in their environment leads to a new equi- 

tory of fundamental experimental and theoretical study, librium molecular weight distribution (MWD). All of this 

is living anionic polymerization jHEllS ,4,5,6, 111111111101 should be compared to conventional inert polymer mate- 

|TT|[T^|T^[ni|T^P^|T7|[TBl[T^|^|niP2H23„24t|^ This rials whose MWDs are frozen, 
class of polymerization is widely employed to manufac- 
ture polymers with nearly monodisperse MWDs and con- 

trolled architectures such as block and star copolymers The subject of this paper is the dynamical response 
Eg. A recent variant on this theme with immense poten- of ^'""'"^S polymers. We study the class of hvmg polymers 
tial technological impact is living free radical polymeriza- ^^^^^ exemplified by hvmg ionic polymerization, having 
tion 27,28 . Other examples include worm-like surfactant following characteristics: (i) monomers add and sub- 
micelles |2Hl3()l;-!ll;-i2l33l:^4l:ihl:ffi| and biological polymers ^ract from chain ends and (u) the total number of chains 
such as actin and tubulin filaments |S71EHlE3EnilIIllll ^^'^ biological world, actm and microtubule 
whose special properties are exploited by living cells for filaments satisfy condition (i), and the dynamics of both 
motility and structural integrity filaments and filament caps are frequently described by 

As polymeric or pofymer-like materials, the unique situation (ii) m Wormlike micelles on the other hand 
feature of these systems is that the chains are dynamic ^O":!*; ^^^ses grow from their ends only |ini|3D| (condi- 
tion (i)) but do not satisfy condition (ii): the number of 

e-mail: bo8@columbia.edu micelles is not conserved. Interestingly, this makes their 

^ e-mail: dv35@columbia.edu dynamical response very different (see discussion section). 
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(a) o + • io 

V. ' 

Fig. 1. (a) Schematic of initiation reaction in ionic polymer- 
ization. Empty (filled) circle represents initiator (monomer), 
(b) Polymerization and depolymerization reactions at the ionic 
chain end of a polymer chain of A*' units. 

Our presentation will focus on living ionic polymer- 
ization, though most of our results are completely gen- 
eral. In these systems, polymers grovif reversibly from their 
ionic ends (see fig. Due to their "living" nature, if 
fresh monomer is added polymers will grow until the extra 
monomer has been consumed; further addition of a dif- 
ferent monomer species, say, will resume polymer growth. 
The main achievement of this technique is a distribution of 
chain lengths with very small Poisson fluctuations around 
the mean. This monodisperse MWD is not however a true 
equilibrium distribution. For sufficiently long times slow 
depolymerization reactions become important and lead to 
a true equilibrium MWD which is in fact extremely broad. 
For many applications, for example styrene, the timescale 
to reach equilibrium is extremely large and depolymeri- 
zation-induced broadening effects are unimportant 4A. 
However, in some cases equilibration times are accessi- 
ble, an example being a-methylstyrene whose equilibrium 
properties have been studied in many experiments 
1^,6,7,8,9,10,11,12,13,. 

The reason that chain length fiuctuations in equilib- 
rium are so large is that even when a state is reached where 
the monomer concentration reaches a steady state value 
such that polymer growth rates exactly balance depoly- 
merization, through polymerization and depolymerization 
reactions monomers continue to be reshuffled among liv- 
ing chains. Thus even if one starts with an essentially 
monodisperse MWD as shown in fig. El^a), random ad- 
dition and subtraction of monomers from chain ends will 
result in a "diffusive" random walk motion of chain ends in 
A'^-space (where N is the number of monomer units added 
to initiator) with a certain characteristic "diffusivity" D. 
Since there is no restoring force to this motion, chain ends 
will eventually diffuse to distances of the same order as 
the mean chain length, resulting in a broad distribution 
as shown in fig.EIb). In a simple mean-field picture where 
monomer-monomer excluded volume interactions are ne- 
glected this leads to the Flory-Schultz equilibrium distri- 
bution gniiiniEi 

0oo(iV) = ^e-^/^- , (1) 

where N^o is the equilibrium mean length. 

The thermodynamic properties of equilibrium living 
polymers have been studied in a large number of experi- 
mental and theoretical works. According to theory, going 




Fig. 2. Fluctuations of living polymer chain length in equilib- 
rium are very large. For clarity, living chains are depicted as if 
extending from an imaginary substrate. Living ends are shown 
as filled circles. Even if one starts with a nearly monodisperse 
MWD as in (a) and with monomer concentration (grey back- 
ground) such that polymerization and depolymerization rates 
balance, random polymerization and depolymerization reac- 
tions eventually broaden the distribution up to distances of 
the same order as the mean A'^, as shown in (b). Note the re- 
sulting MWD includes chains of "zero" length (free initiators). 



beyond mean field by accounting for excluded- volume in- 
teractions leads to power law modifications jl8lll9lE(11l34| 
to eq. Moreover, an analogy jl4lll5l[TI)lll7| between a 
living polymer system and the much studied Ising spin sys- 
tem in the formal limit of spin dimensionality n becoming 
zero (similar to the well-known analogy in the thermody- 
namic theory of semidilute polymer solutions |47l48j which 
can in fact be formulated equivalently to a living polymer 
system showed that the polymerization phase transi- 
tion near the polymerization temperature is second-order 
with a number of characteristic power-law exponents. The 
results of a large number of experiments by the group of 
Greer, measuring both the MWD itself 13 and other ther- 
modynamic properties 2 4 are apparently consistent with 
the above theories, favoring somewhat the non mean-field 
approach. 

The most essential feature of living polymers however, 
namely their ability to respond dynamically to external 
perturbations, remains largely unexplored. How does an 
equilibrium living polymer system respond to a small per- 
turbation? A conceptually simple perturbation is the addi- 
tion of a small amount of extra monomer (an "m-jump"). 
A similar perturbation, easier to perform in practice, is a 
"T-jump," i. e. a small sudden change in temperature. Un- 
der such a perturbation, how rapidly and in what manner 
will the monomer concentration, MWD and mean chain 
length reach a new equilibrium? In this paper we address 
these questions theoretically. We will compare our pre- 
dictions in section 8 with small T-jump viscosity relax- 
ation experiments by Greer et al. |11) . This same group 
has also monitored monomer and MWD dynamics fol- 
lowing much stronger perturbations, finite temperature 
quenches [T^IOj (see fig.OJ. 

Our main broad conclusion is that, counter- intuitively, 
the response to a very small perturbation is extremely 
strong: even though the initial and final MWDs (the latter 
attained after relaxation is complete) are very close to one 
another, the MWD at intermediate times is very strongly 
perturbed from equilibrium, i. e. its shape is very different 
from either the initial or final equilibrium shapes. In this 
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Fig. 3. Schematic of one of the classes of experiment performed 
by Greer's group [T^IT^ . The mix of monomer, initiator and 
solvent is stored for sufficiently long above the polymerization 
transition temperature Tp such that all initiator is activated. At 
this point all living chains are of "zero" length; we refer to these 
as "free initiators" (bifunctional dimers in this case). After 
a quench below Tp, polymerization onsets and the system's 
response to abrupt temperature changes is then studied. In 
this paper we study more elementary perturbations: very small 
T-jumps of an already equilibrated system. 



sense, there are no small perturbations: a small thermo- 
dynamic perturbation well-described by a linear suscepti- 
bility (i.e. one inducing a small change in the equilibrium 
state) has nonetheless a large dynamical effect; intermedi- 
ate states deviate strongly from equilibrium in the sense 
that observables are perturbed in a highly non-linear man- 
ner. We refer to this as dynamical ultrasensitivity. 

A hint of this sensitivity is already apparent in equilib- 
rium, where a very delicate balance is established between 
growth and shrinkage such that the net polymerization ve- 
locity V is essentially zero, of order 1/ N^o (due to the small 
fraction of special non-depolymerizable chains of length 
TV = 0). Now the effect of a small T-jump (say) is to pro- 
duce a small but finite velocity which overwhelms this tiny 
velocity value and destroys the delicate balance sustain- 
ing equilibrium. Thus if for example v becomes positive, 
the MWD will uniformly translate toward larger molecular 
weights, creating a depletion in the MWD at small chain 
lengths, as shown schematically in fig0| This translation 
process stops only after sufficient monomer has been poly- 
merized such that the monomer concentration drops close 
to its destination equilibrium value, corresponding to a re- 
covery of the zero net growth rate equilibrium condition. 
We find that only on a much longer timescale does "dif- 
fusion" of chain ends fill up the hole and lead the system 
to the destination equilibrium MWD, as shown in fig. 01 

Why does this hole develop in the MWD? The rea- 
son is that chain length scales homogenized by diffusion, 
(-Dt)"'^/^, are smaller than scales affected by translational 
motion, vt, for times longer than the timescale t* « Djv^ . 
Thus if t* is shorter than the timescale needed for transla- 
tion to stop, a depletion hole develops in the MWD. But 
even when translation is completed earlier than t*, we 
find that the perturbation is still nonlinear; even though 
no hole is formed there is still a large amplitude reduction 
in the MWD, much larger than the relative magnitude of 
the initial perturbation, e. 



This work is related to earlier theory of Miyake and 
Stockmayer |44j who, following a prior treatment by Brown 
and Szwarc |50| . studied analytically and numerically the 
dynamics in the special case where initially all monomer 
is unpolymerized (see fig. [SJ. Their analytical results, ap- 
plicable in the case of very small depolymerization rates, 
showed a 2-stage relaxation. During the first, lasting a 
time ~ iVoo , living chains grow coherently from zero length 
up to the equilibrium value. In the second stage the shape 
of the sharply peaked MWD relaxes slowly toward the 
equilibrium exponential distribution. Although their per- 
turbation analysis broke down during this stage, Miyake 
and Stockmayer estimated that its duration scales as . 
Later theoretical studies addressed the second stage dy- 
namics 51,52, 53 of this special case. 

A short announcement of the present work has al- 
ready appeared jSl]. The structure of the paper is as 
follows. In section 2 we establish the dynamical equa- 
tions obeyed by living polymers which show the nonlinear 
coupling between the monomer population and the living 
chain MWD. For simplicity, we neglect excluded volume 
interactions and chain length dependent polymerization 
rates. Starting from the equations of section 2 we then 
analyze the response to a T-jump. This response depends 
on the sign of the velocity induced by the perturbation. 
The case of positive initial velocity is analyzed in sections 
3 and 4, while the negative in section 5. The special case 
of very small T-jumps is analyzed in section 6. In section 
7 we show that the results of the previous sections di- 
rectly generalize to small perturbations of arbitrary form. 
Finally we conclude with a discussion of the results, and 
the experimental outlook. 

2 Living Polymer Dynamics: 
Monomer-Polymer Coupling 

An important aspect of living polymers dynamics is the 
coupling between the two distinct species present in the 
solution, free monomers and living chain MWD: on the 
one hand living chain growth rates depend on the concen- 
tration of free monomer and on the other hand, due to 
mass conservation, monomer concentration if a function 
of the living chain MWD. Stated differently, living chains 
grow according to a velocity field which is self-consistently 
updated depending on the response of the MWD. In this 
section we develop the equations obeyed by monomer and 
MWD. The dynamical response to perturbations is dis- 
cussed in the following sections. 

A crucial feature of the living polymer systems we 
study is that the number of living chains remains fixed. 
This is what arises in the experiments of Greer et al. 
OSEli where this number is determined by the amount 
of initiator initially added in the solution. Denoting zq the 
concentration of living chains (i. e. chains which have been 
initiated), mass conservation implies that monomer con- 
centration, mj, and MWD, (f>t{N), [t denotes time) obey 

/•oo 

mt = mtot - i^Nt , Nt= dN N (pt{N) , (2) 

"'0 
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Fig. 4. There are three stages in the relaxation of living polymers to equilibrium. During the first, living polymers grow 
coherently as in (a). The translation stops when they have consumed enough monomer such that monomer concentration 
drops to a value for which polymerization and depolymerization balance (shown in (b)). During the second stage, (c), through 
"diffusive" random polymerization and depolymerization reactions, the MWD slowly fills the depletion of chains of short lengths. 
The third and final stage entails global shape relaxation of the MWD on a difi'usion timescale corresponding to Noo ■ 



where mtot is the total monomer concentration (includ- 
ing polymerized monomers), and Nt is mean chain length. 
Here (j)t is normalized to unity and N — corresponds to 
free initiator. The coupling between monomer concentra- 
tion and the MWD is manifest in eq. 

Now the dynamical equations equations obeyed by (j)t 

are 



dMN) _ djt{N) 



dt 
Jt{N) 



dN 

VtMN) - 



ON 



Jt(0) = 0, (3) 



where 



rrit 



{k+mt + v-)/2 . (4) 



Here fc+ is the propagation rate constant and is the 
depolymerization rate. The "diffusion" coefficient Dt, de- 
scribes fluctuations in the rate of polymerization/depoly- 
merization. Eq. (O is identical to the diffusion dynam- 
ics of a particle in a linear potential with time-dependent 
slope —ft whose time-dependent diffusion coefficient is Dt 
(see fig. El . The novelty here is that due to the coupling 
between the MWD and monomer, both the slope and dif- 
fusivity are functionals of (j) itself. 

The reflecting boundary condition in eq. © at TV = 
represents the fact that when a living chain looses all 
of its monomers and becomes a free initiator molecule 
(i. e. reaches length A'^ = 0) it cannot depolymerize fur- 
ther and must grow again. Now we restate eq. lO as 



dt 



-vt 



dN 



D.. 



92 



0t(O)M(O) = AM , (5) 



where = d<f)t/dN. 

Now it might at first seem that despite the coupling 
of eq. (O, the monomer subsystem would effectively un- 
couple from the MWD dynamics, since one might naively 
expect that the monomer dynamics do not depend on 
the shape of the polymer MWD but only on the num- 
ber of living chains. This would imply simple first order 



kinetics for m* resulting in exponential relaxation of the 
monomer concentration, independent of the MWD dy- 
namics. The coupling however arises from the existence 
of special chains in the MWD of zero length, i. e. free ini- 
tiator molecules of which there are (^t(O), which unlike all 
other chains cannot depolymerize. Indeed, the dynamics 
obeyed by mt, or equivalently Vt, are 



dvt 



Tfast 



A0t(O) 

Tfast 



Tfast = 



/ 



(6) 

We see clearly that the only aspect of the living MWD 
which monomers see is 0t(O). Eq. (jSJ is derived by calcu- 
lating dNt/dt by multiplying eq. © by N, integrating over 
all N , using the reflecting boundary condition, and then 
using eq. The relationship between Tfast and N^o fol- 
lows from eq. (jHJ below. Here we have introduced the two 
basic independent dimensionless parameters of the sys- 
tem: 

f = TT , r = . (7) 

k+mtot mtot 

The physical meaning of / is the following: if the sys- 
tem were a pure solution of unpolymerized monomer (plus 
solvent and initiator) / would be the ratio of backward 
to forward polymerization velocities. The value of / is 
temperature dependent, being unity at the polymeriza- 
tion temperature and smaller or larger than unity in the 
polymer and non-polymer phase, respectively. The param- 
eter r, namely the ratio of living chain to total monomer 
concentration, is independent of temperature and is al- 
ways much smaller than unity. Its smallness is related to 
the mean degree of polymerization being much larger than 
unity in equilibrium (see eq. © below). 

Now setting the time derivative in eq. JSJ to zero and 
using eq. (0) it is shown in Appendix A that in equilibrium 
the Flory-Schultz distribution of eq. is recovered and 
one has to leading order in r: 



l-f 



mtot 



f 



V 



D^ 



(8) 
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Fig. 5. Schematic of tiie MWD dynamics as analyzed by 
Miyake and Stockmayer |44f and their theoretical predictions 
regarding the timescales involved, (a) Initially {t — 0) all ini- 
tiators are unpolymerized. (b) For t <^ rfast chains grow coher- 
ently with a narrow MWD until its peak reaches the equilib- 
rium mean length Naa at times of order Tfast • (c) During a third 
stage lasting up to Taiow, which was beyond their theoretical 
analysis, the MWD spreads to the equilibrium MWD. 



where r ^ 1 and subscript oo denotes the t ~* oo equi- 
librium value for the corresponding variable. We assumed 
(1 — f)lr^^'^ ^ 1, i.e. that the temperature is not ex- 
tremely close to the polymerization temperature. 

An important feature in eq. JHl is that in equilibrium 
the velocity, i. e. the slope of the "potential" term in eq. 
JSJ, settles down to a very small negative value. Were there 
no diffusion, living chains subject to a negative velocity 
field would shrink to zero length. However due to the small 
magnitude of the velocity, diffusion is strong enough to 
broaden the MWD up to distances of order Nqo-, as illus- 
trated in fig. 121 We show in the following section that any 
apparently small external perturbation has an enormous 
effect on the velocity which, depending on the direction 
of the perturbation, may become either very negative or 
very positive. 



3 Response to T-jump; Positive Velocity 
Boost 

Perhaps the simplest way to perturb an equilibrium liv- 
ing polymer system is by a small temperature change (T- 
jump). Taking as an example a-methylstyrene, the data of 
fig. show equilibrium fraction of polymerized monomer 
as a function of temperature; evidently, changes in T lead 
to different values of the equilibrium monomer concentra- 
tion moo and therefore of as well. In this section we 
consider relaxation dynamics after an equilibrated system 
is subjected at t = to a small temperature change 8T 
such that the system will eventually reach new equilibrium 
values TOoo and TVqo. That is, we follow the dynamics of 
the transition from an old equilibrium towards a slightly 
different new equilibrium state. 

The magnitude of the perturbation is measured by the 
small parameter 

6^^, (9) 



^(N) 




Fig. 6. Living polymer dynamics are equivalent to diffusion in 
a potential with slope —vt- In equilibrium the slope reaches a 
very small value of order 1/iVoo (solid line) resulting in a broad 
MWD. Under a small T-jump the slope is strongly perturbed 
and depending on the direction of the jump becomes either 
positive (dashed line) or negative (dotted line). 



The value of e is related to the magnitude of the T- 
jump. For example in the case of a-methylstyrene, assum- 
ing that the system is initially below the polymerization 
temperature Tp in the region of fig. [3 between 260° K and 
280°K, one has using eq. Q 



(5T 
'50°K 



(a-methylstyrene) 



(10) 



Thus for a-methylstyrene a temperature increase, bT > 0, 
results in a negative ^ttiq, i.e. a reduced initial monomer 
concentration relative to the destination equilibrium value. 
Now the i = relative perturbation of the MWD is 



This is of order e and smaller than unity for all N 
from the unimportant values N 3> Nqo where the 
is exponentially small) as shown in fig. IHl^a). 

Because velocity, monomer concentration and 
chain length are all linearly related (see eqs. (01, 1 
Q) changes in these quantities are simply related 



1-/ 
/ 

(11) 
(apart 
MWD 



mean 
1) and 



as 



Smt 



SNt 



Svt 



(12) 



These relations are true for all times. They allow us to 
follow the dynamics of velocity alone. Once vt is known, 
Nf, mt, Dt are determined. 

From eqs. H12|l and the initial relative changes in 
V, D are: 



5vq 



SDq 



(13) 



where we define 5mt = rrit — rrioa and similarly for other 
quantities. All equilibrium values refer to the destination 
{t — oo) equilibrium values. Thus Smo is the initial (t = 0) 
deviation from the final equilibrium. 



The important point is that since iVoo ^ 1, the relative 
change in v is much larger than e : the velocity is highly 
perturbed, as shown in fig. El Depending on the sign of e, 
the velocity may remain negative as in equilibrium, or it 
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Fig. 8. (a) Sketch of the initial and final MWDs during a T- 
jump experiment. Both distributions are exponential. Though 
very close to one another, in the course of relaxation the time- 
dependent MWD becomes very different to the final exponen- 
tial MWD. (b) General perturbation (non-exponential initial 
MWD). Most T-jump results derived here apply to this general 
case also. 
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Fig. 7. Fraction of polymerized monomer in equilibrium as a 
function of temperature for living a-methylstyrene in tetrahy- 
drofuran solvent, initiated by sodium napthalene. Data from 
ref. Ratio of mole fraction of total monomer, mtot, to ini- 
tiator is 0.0044. Mole fraction of total monomer is 0.1538. The 
polymerization temperature Tp is near 295°K. 



may change sign. Since the relative perturbation in D is by 
contrast small, we see that the delicate velocity-diffusion 
balance which sustained equilibrium is now destroyed. 

We will show below that relaxation to equilibrium now 
occurs in three stages. During the first, dominated by 
translational motion, the MWD is boosted far from its 
equilibrium shape. The next stage involves diffusive restora- 
tion of the region of the MWD which suffered maximum 
distortion during the first stage. The third and final episode 
entails a very slow diffusion-dominated recovery of the 
global MWD. For the remainder of this section we treat 
the case of a positive initial velocity boost, where chains 
initially grow {6T < for a-methylstyrene). Negative ve- 
locity boosts, where chains initially shrink, are considered 
in section 5. 



3.1 Coherent Chain Growth: Hole Development 

During the first stage of the relaxation process, velocity 
dominates over diffusion since v has been so strongly per- 
turbed. Chains thus grow coherently, consuming the ex- 
cess monomer in a timescale rfast by which time transla- 
tion will have halted and a highly non-linear hole will have 
developed in the MWD (see fig.EJ. Here Tfast, defined in 
eq. © , is the time for the MWD to translate distance SNq 
and reach the destination mean, i.e. Tfast ~ ^Nq/Svo (see 
eqs. CSIl and using eq. (jHJ). 

To see all of this more quantitatively, consider the ve- 
locity dynamics eq. JBJ. Initially, the velocity term on the 
rhs is much larger in magnitude than the diffusion term. 
We show in Appendix B that this remains true up until a 
time Tqs, defined below. It follows that for these times Vt 



relaxes exponentially. The same is true of mt,Nt and Dt 
which we recall are linearly coupled to Vt (see eq. ((T^ ). 
Thus 



Svt 



ev e 



(14) 

Note that for t ^ Tfast there are no free initiators (at 
much longer times these are restored, see below). Thus 
for these intermediate times all chains are identical as far 
as on and off monomer kinetics are concerned, and the net 
rate of monomer-polymer mass exchange can only vanish 
if Vt vanishes. For this reason vt decays to zero, though 
much later it will recover its small negative equilibrium 
value. The 4>t dynamics, eq. (jSJ, thus simplify to 



dt 



o(Pt 
ON 



(i«Tqs) (15) 



where the diffusion coefficient may be approximated by 
Doo since the contribution of 6Dt is negligible (see Ap- 
pendix B). This linear equation, plus the time-dependent 
boundary condition of eq. Q , has solution 



dN'G'r''{N,N')MN') (i«Tqs) (16) 



which describes a translating and simultaneously broad- 
ening MWD, as shown in fig.El Here Gf'''''{N,N') is the 
propagator of eq. H15() whose properties are calculated in 
Appendix B. 

An important quantity in what follows is (/'t(O), i. e. the 
concentration of free initiators. Using eq. p6|l we show in 
Appendix B that 



0t(O)/</)oo(O) « 

F{t*/tfl^e-*/'' {t* «i«Tfast) 



4D„ 



(17) 



where C = 4/7r^/^ and F ^ tt"^/^. Eq. I^Tj) has a clear 
physical meaning: as the MWD translates, the position 
of its peak at A'' = after time t has shifted to chain 
lengths of order vot. Were diffusion absent, this would have 
created a depletion region in the MWD at small chain 




N„ Noo No Noo No Nco 

(a) t«t* (b) f*« t « Tfo,t (c) ^ = 'Cfast 



Fig. 9. Initial response of tlie MWD after a positive boost, case of strong perturbation (t* <^ ifast)- (a) For small times, chain 
end "diffusion" smooths out the hole which would otherwise form at small due to chain growth, (b) At larger times coherent 
chain growth beats diffusion and a hole develops in the MWD. (c) Coherent growth stops when mean chain length reaches a 
value very close to Noo, displaced by an amount eNoo/d from the original value TVq. For typical cases, the parameter 6 is of 
order unity. 



lengths. However diffusion smooths out inhomogeneities 
on distances of size {Doct)^^^- Hence for times shorter 
than the crossover time t* for which vot* ~ {Doot*y^^ , 
diffusion has enough time to fill the translationally in- 
duced hole and the relative deviation from equilibrium is 
small. For times longer than t* however, the MWD trans- 
lates a distance away from the origin much larger than 
what diffusion could have homogenized, a hole develops 
in the MWD, and the concentration of zero length chains 
becomes exponentially small. 

In eq. H17|) we assumed that the perturbation is large 
enough such that t* <C Tfast • The special case of extremely 
small perturbations where this is no longer true is dis- 
cussed separately in section 5. 

3.2 Diffusive Length Relaxation: Hole Filling 

For times t 3> Tfast the fast variables monomer and mean 
chain length have relaxed very close to their equilibrium 
values. We have seen that velocity becomes exponentially 
small. Meanwhile (see below) the number of free initiators 
4>tiO) is gradually recovering. At a certain moment, there- 
fore, the 2 terms on rhs of eq. © become equal to one an- 
other and we show in Appendix B that for all later times 
the velocity time derivative is much smaller. Hence the 
velocity dynamics are now fundamentally changed. The 
new regime is one of quasi-static evolution, enslaved to 
the dynamics of 0t(O): 

Svt^ -V'S(f>t{0) (t>Tqs) . (18) 

Physically, this reflects the fact that the the only aspect of 
the MWD shape seen by monomers is the amount of spe- 
cial non-depolymerizing zero-length chains. We see that 
as the MWD slowly rearranges itself, so the fast variables 
vt , rrit and Nt variables respond quasi-statically. 

Notice that the quasi-static regime does not onset af- 
ter Tfast but rather after Tqs. The reason is that during the 
initial boost the term DtiptiO) in eq. © decays exponen- 
tially on a timescale t*; its magnitude therefore at Tfast 
is much smaller than vt which decays on a timescale Tfast 



(see eqs. H14|) and (I17|| V There is thus a cross-over period 
before the quasi-static balance between the two terms is 
established. 

Now let us examine the MWD dynamics. Since vt be- 
comes exponentially small after the MWD stops translat- 
ing, we have 

d(hf d'^chf , X , s 

^ « i^oo ^ (Tfast « t « Tfill) (19) 

with reflecting boundary conditions at the origin. Here 
once again the contribution of &Dt is negligible (see eqs. 
(I12|l . 118|l . and (|21ll ) and the corresponding term has been 
neglected in eq. (|19|) . Here 

Tfill = ^ Tslow , Tslow = . ^ ■ (20) 

The timescale TsIow is the longest relaxation time of the 
system, the time for the slow global shape characteristics 
of the MWD to relax. It equals the diffusion time for the 
MWD width Noo- Meanwhile Tfm is the diffusion time for 
the hole width, eNoo/0- It is shown in Appendix B that 
eq. H19|l has solution 

1 /2 

^ H e-^«"/* (Tfast « t « Tfiii) (21) 

where H — vr^^/^. Thus for times shorter than Tfin the 
concentration of zero length chains remains exponentially 
small while for longer times the hole fills (see fig. I10|l and 
(/)f(0) recovers its equilibrium value. 



3.3 Linearized Dynamics at Long Times 

For times beyond Tfin the hole-filling is complete and the 
MWD's non- linear feature has disappeared. Thus, finally, 
a truly linear regime onsets: relative deviations of all vari- 
ables from equilibrium are less than unity and perturba- 
tion theory can now be applied. This is done in Appendix 
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^ « -SMO) « 

V 

J0oo(O) (TfinA)l/2 

Me0oo(O)(t/T,iow)- 



3/2g-t/(16T,io 



lo 

{t > "^slow) 



(22) 



where J, M are positive constants of order unity. Thus 
(50t(O), and all of the fast variables which remain quasi- 
statically enslaved to 0t(O), relax exponentially on a time- 
scale Tsiow with a power law prefactor. 



4 The Second Moment of the MWD Relaxes 
More Slowly than the First 

A picture has emerged of fast mass- related variables which 
relax in Ttast and slow MWD shape properties which relax 
in the much longer time Tgiow The MWD's first moment 
defines the total mass in the polymer system: this is a fast 
variable and was considered in section 3. The simplest slow 
shape property is the dispersion At, 



At EE Ni 



dNN^MN) , (23) 



which is closely related to the second moment. In this 
section we consider the very different and slow relaxation 
of this quantity. 

Substituting eq. Hll|) in eq. (|23|l one finds that fol- 
lowing a T-jump the initial relative perturbation in the 
dispersion, SAq, is 



^^0 _ g £ 
A^ 9 ' 



(24) 



where the equilibrium value, Z\oo, is found using eq. (|T)). 
Now the dynamics of the dispersion can be derived by 
evaluating the time derivatives of N^, Nt in eq. (|^ . These 
can in turn be calculated by multiplying eq. ^ by N'^ 
and N, respectively, integrating over all N and using the 
boundary condition of eq. JSJ- One has 



1 d6At 
dt 



2r. 



low 



Do 



Nt 0t(O) 

A^oo0oo(O) 



1 ^0t(O) 



2tsi, 



ow V^ctt 



(0) 



(25) 



Here we replaced Dt and Nt with their t = oo values. This 
is correct to leading order since the relative perturbations 
in these quantities are much smaller than that of '?!>t(0) 
which undergoes much larger changes (see section 3). 

Now integrating eq. 1(23) using eqs. (ETj), ijllj and if^ 
one finds that up to the time Tf^\\ the dispersion has changed 
very little relative to its value just after the perturbation 
at t = 0: 



At-Ap 
A^, 



r(C/3)(t/r3iow)(tA*) 

\t/(2Tslow) 



1/2 



{t < t*) 

{t* -^t-t: Tfill) 



(26) 



^(N) 



m- ( 




Fig. 10. Schematic of MWD during hole-filling episode. Dot- 
ted (dashed) lines show MWD profile for times rfast ^ i ^ Tmi 

(Tfill < t < Tslow). 



This is in stark contrast to the first moment which had 
already relaxed by the much shorter time Tfast- The dis- 
persion's relaxation process does not properly get going 
until times of order TsIow: 

SA ^ ( SAo/A^ + Je6'-i(t/T,iow)^/^ (Tfiii < t < t,iow) 
A^ ~ \ae (VT-siow)-'/'e-*/(if^"='-) [t » t^iow) 

(27) 

where a = 8A//39. The relative perturbation in the dis- 
persion remains of order e for t <C Tgiow and subsequently 
relaxes exponentially. 



5 Response to T-jump; Negative Velocity 
Boost 

In this section we study the case opposite to section 3, 
when e < and the T-jump induces a higher equilibrium 
monomer concentration and a smaller N^c (see eq. ©)• 
Such a perturbation in the example of a-methylstyrene 
would be produced by an increase in temperature as one 
may see in fig. [7| In this case the sign of the initial per- 
turbed velocity (eq. (|13|l ) is negative (see fig. (HI). Thus 
initially the MWD is not boosted towards larger molecu- 
lar weights as in section 3, but is instead boosted in the 
opposite direction, towards smaller N. Therefore, instead 
of depletion, an excess of small length chains and free ini- 
tiators is produced, creating instead of a hole, a sharp 
peak in the MWD near = as shown in fig. ^2 

Despite the differences between the two cases, many 
of the results of section 3 remain unchanged. In fact all 
equations of section 3 up to eq. (|16|l remain unchanged, 
the only difference being that e is negative. Thus eq. (|16|l 
now describes a MWD coherently shrinking and simul- 
taneously broadening with diffusion coefficient Doo- But 
since free initiators cannot depolymerize, excess polymer 
must build up near N = 0. Indeed, starting from eq. 
it is shown in Appendix D that 



>o(0) 



F' [t/t*) {t* <^t^ Tfast) 



(28) 
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Noo Nn 




No, No 




No, Nn 



(a) t «t' 



(b) r«i«Tj 



fast 



fast 



Fig. 11. Initial response of MWD after a negative boost (strong perturbation, t* <C Tfast)- (a) Chains shrink coherently and 
MWD moves from the initial mean iVo towards iVoo- For small times, diffusion smooths out excess peak at small A'^. (b) Beyond 
t* diffusion is too slow to counteract coherent chain shrinkage and short chains accumulate at the origin, (c) Coherent shrinkage 
halts when mean chain length approaches Noo- 



where C" = 4/71^/^ and F' = 4. Eq. has a similar in- 
terpretation to eq. H17|l . For t <^t* diffusion is fast enough 
to smooth out the excess polymer accumulated by the neg- 
ative velocity at the origin and thus 0t(O) remains close 
to its initial value (see fig. ^2 (a))- For t ^ t* however the 
flow towards the origin is faster than what diffusion can 
smear out and a peak forms whose height increases with 
time, as shown in fig. ITIT b). (Once again it is assumed 
that the perturbation is not so small that t* <C Tfast; the 
opposite case is examined in section 5). 

After the completion of the first stage at times of order 
Tfast) similarly to section 3, the fast m^, Nt variables have 
essentially relaxed and start to respond quasi-statically 
to the slow diffusion-driven shape changes of the MWD. 
In fact we show in Appendix D that eqs. H14() and H18|l 
still apply, but now with a different crossover timescale, 
Tqs = Tfast ln ^oo ■ This difference arises because unlike the 
positive boost case (/'t(O) does not become exponentially 
small and so needs less time to catch up with the vt term 
in the velocity dynamics, eq. JBJ. 

Similarly to the positive boost, eq. (|19|l still applies 
for t < Tfiii. By this time the accumulation at the origin 
is able to diffuse to distances of the same size as the re- 
gion from which it was transferred from during the first 
stage. Hence (f>t{0) becomes of order its initial value, and 
thereafter perturbation theory is valid. In Appendix D we 
show that 



6t(O)/0oo(O) « 

1 + J'(Tfiu/t)i/2 
H-M'e(t/Tsiow) 



-3/2 -t/(16r,io„) 



(Tfast < t <C Tfiii) 
(Tfill < t < Tslow) 
{t > Tslow) 

(29) 



where K',J', and M' are positive constants of order 
unity. Note there is apparently no smooth cross-over at 
t = Tfast between the forms of eqs. (|28|l and H29|l. This 



is because there is in fact an extra episode at times of 
order Tfast during which the very thin peak spreads from 
width {Doot*Y^'^ to a width (ZJcxjTfast)^/^- Note also that 
the time behavior for t 3> Tfm has been derived using the 
linearized dynamics of Appendix C. 

Finally, to obtain the dispersion dynamics we use eqs. 
and (|2ni) in eq. to obtain 



SAt/A^ - SAo/Aoo ~ 

-(C73)(t/Tsiow)(t/t*)i/' 
-F74(t/Tsiow)(Vi*) 

-SAo/A^ - (8Af 739) e (ts. 



{t < t*) 

(i* < i < Tfast) 
(Tfast < t < Tfill) 
(Tfill < t < Tslow) 

3w/t)^/2e-*/(i6-.iow) 

(t » Tslow) 



(30) 



which is very similar to eqs. (|26|) and (|27|l for the positive 
boost case. Again the dispersion relaxes on a timescale 
Tslow , much longer then the relaxation time of Nt . 



6 Very Small T-jumps: e less than its critical 
value ec 

Hitherto we implicitly assumed sufficiently strong velocity 
boosts that the timescale t* after which translational mo- 
tion dominates diffusion and a hole or peak begins to form 
is shorter than the time for full development of the hole 
or peak, namely Tfast • Now the ratio of these timescales is 



t* 



Tfast 



9I/2 



1/2 



(31) 



which defines a critical value of the perturbation param- 
eter, Ec- Clearly, the system's response must be different 
for e values below this rather small threshold. We consider 
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f(0) 



\\ 80(0) - (^) ' 




Fig. 12. Schematic of MWD development during rfast <C t >C 
Tsiow when the small chain peak spreads diffusively. Dotted 
(dashed) lines correspond to times rfast <^t<^ rmi {Tmi -^t <^ 

rslow ) • 

such very small perturbations in this section. Although a 
full hole or peak does not have time to develop, we will 
find that the response is still strong and non-linear. 

Since such weak perturbations do not involve the cre- 
ation of a complete hole or a peak in the MWD, one ex- 
pects that the timescales t* and Tm associated with the 
their creation and destruction lose their physical mean- 
ing as cross-over timescales. This is shown in Appendix 
E where, using similar arguments to the ones of sections 
3 and 5, it is proved that weak perturbations are identi- 
cal to stronger perturbations except that the regime up 
to t* is truncated at rfast and the following regime (which 
for stronger perturbations lasted till rfin) deleted. In Ap- 
pendix E it is shown that for both positive and negative 
boosts, one has for the fast variables 



Svt 



e-tM„3t (t< rfast) 
-'5<?!>t(0) (rfast <i<Ts 

low J 



(32) 



Note that that the timescale rqs does not appear since its 
magnitude lies between those of t* and Tfm. 

The evolution of the MWD (see figs.[l21lllll is similarly 
shown in Appendix E to lead to the following solution for 
0f (0) in the case of a positive initial velocity: 



l-C{t/t*y/^ (t«Tfast) 
1 /2 

1 - Q (Tfill/i) (rfast < i < Tslow) 



{vo > 0) 



while for a negative initial velocity one has instead 

hiO) f l + C'(t/r)l/2 (t« rfast) 



(33) 



1 + Q'(Tfin/t)i/2 (rfast < t < rsiow) 



{vo < 0) 

where Q, Q' are positive constants. Finally, it is straight- 
forward to show that the relaxation for for t 3> rsiow is 
exponential as for the stronger perturbations, eq. H22|l . 



The most important result of this section is that even 
though the e < ec case is "weak" as compared to sections 3 
and 4, the system's response remains large and nonlinear. 
One sees from eqs. (T^ and (p^ that at rfast when the 
deviation of 4>t{0) from equilibrium is largest, the relative 
deviation of 0t(O) from equilibrium is much larger than e. 
Therefore even in the weak case a hole or peak does form 
in the MWD but of a smaller amplitude as compared to 
the stronger perturbation case (see figs. lT^IT^ . 

Now the dispersion dynamics may be similarly evalu- 
ated using eqs. (|33|) and (|34|l in eq. (|25|l . One finds that 
At relaxes once again on a timescale rgiow ■ For the positive 
initial velocity case one finds 



SAt/Aoo - SAo/Aoo « 

r (C/3)(t/rsiow)(tA*)i/' 
I Oe0-i(</rsiow)i/2 
I -5zio/Z\oo + (8Af/39)e 



[t < rfast) 

(rfast <^t^ rsiow) 

rsiow/t)^/2e-*/(i6-=- 



(35) 



(t 3> rsiow) 
while for a negative initial boost 



\l/2 



SAt/A^ - SAq/A 

-(C'73)(t/rsiow)(t/i 
-Q'e0-Ht/rsiow)i/' 
-SAo/A^ - (8AfV39) e (rsiow 



(t < rfast) 

(rfiii < t < rsiow) 



/t) 



5/2, 



{t > rsio 



) 



7 General Perturbations 

We have so far considered a specific type of perturbation, 
namely T-jumps. In this section we consider general per- 
turbations. We will see that almost all of the phenomenol- 
ogy remains unchanged. 

Let us start with another simple type of perturba- 
tion: addition of a very small amount of monomer to an 
equilibrated system. (This would require special care to 
avoid introducing impurities which in current experimen- 
tal methods are destroyed by the initiators before poly- 
merization |2|.) Changing the total amount of monomer 
perturbs both independent system parameters, r and / 
(see eq. JTl), while a T-jump perturbed only /. Despite 
this, the results of the previous sections remain identical. 
The reason is that since in both cases the initial MWD is 
an exponential distribution whose mean is different from 
the new target equilibrium MWD, the only variable pa- 
rameterizing the perturbation is SNq, or equivalently SrriQ. 
Thus one sees from eq. Q that monomer addition is sim- 
ply equivalent to a positive velocity boost induced by a 
T-jump. 

Consider now the most general perturbation, i. e. one 
generating an initial MWD of arbitrary shape as shown in 
fig.|Sfb). (This should be compared to the T-jump and m- 
jump which produce initial MWDs of exponential shape 
with a mean slightly displaced from equilibrium.) Such a 
general perturbation is small provided 



\SMN)\ 

0oo(^) 



(37) 
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N„ N„ 



(a) t«Xf., 



N„ N„ 
(b) t = X,i 



Fig. 13. Weak perturbations, e < e^. For such small perturba- 
tions t* ^ Tfast- (a) The coherent component of chain growth 
shifts the MWD mean towards Noo- Diffusion is dominant, 
smoothing depletion at small A^. The net effect is a nonlin- 
ear deviation from equilibrium, (b) By rfast coherent growth 
halts, when mean chain length approaches N^o- The behavior 
for longer times, t S> Tfast, is essentially the same as that shown 
by the dashed lines in fig. ^| 



for all N . As for the T-jump case, we have now 



5mQ 



5No^ 



N, 



oo Jo 



dNNSMN) (38) 



Since eq. 1)38(1 is of the same form as eq. all the argu- 
ments of section 3 leading to the large relative perturba- 
tion in V (eq. ((13|l ) remain unchanged. Moreover, in our 
analysis of sections 3, 4, and 5, the exact shape of the per- 
turbed MWD did not matter in any of the leading order 
results we obtained. The reason is that the effect of the 
velocity field on the MWD is so drastic that to first order 
the response of any MWD close enough to equilibrium is 
the same, in the sense that the MWD simply translates, 
creating a hole (or peak). Thus the analysis of sections 3, 
4 and 5 directly apply to the general case. 

A special exception is when the first moment of the 
perturbation S(f)t{0) is arranged to have a value much 
smaller than eNoo- In this case the relative amount of 
monomer-polymer mass transfer is much less than e and 
this in turn reduces the velocity perturbation. For exam- 
ple, if the first moment were chosen to be so small that 
^mo/moo ~ l/A"oo, then from eq. H13|l one sees that the 
relative perturbation in vt is of order e, much less than the 
order eiVoo produced by a T-jump. In such cases, linear re- 
sponse applies for all times. However, if the first moment 
of S(f>o is much greater than unity the response remains 
nonlinear. 



8 Discussion 

The peculiarity of living polymers is that different mo- 
ments of the MWD relax to equilibrium on different time- 
scales. In the present work we have shown that the first 
moment, i.e. the mean length, Nt, relaxes on a timescale 
Tfast- However the shape of the MWD measured by the 
second moment, or equivalently the dispersion. At, relaxes 
after TsIow Because the two timescales depend on different 



,1 5(|)(0)~(f/f*j" 




N„ No 



(a) t«Xf, 



N„ No 
(b) t = Xf,„ 



Fig. 14. As in fig. ITsl but now for a negative initial velocity 
boost; a peak is formed instead of depletion. Evolution for 
t 2> Tfast as dashed lines of fig. fT^ 



powers of Noc their ratio may be extremely large: 



Tfast 



/ 



V- (l-f) ' 



Tslo 



iV2 

4v- 



(39) 



where l/v~ is the natural time scale and / « "^oo/"^tot is 
approximately the fraction of non-polymerized monomer. 
Our only assumption has been that / is not extremely 
close to unity, i. e. the system is not very close to the 
polymerization temperature. 

In addition to the MWD, the other important observ- 
able is the free monomer concentration, mt- We emphasize 
that since the number of chains is fixed by the number of 
initiators, thus mt is linearly coupled to Nt by mass con- 
servation for all times. Hence the relaxation of these two 
quantities is identical and mt relaxes after Tfast • Note that 
for certain other living polymers such as end-polymerizing 
wormlike surfactant micelles, this is not true |29I3(JI5!3] . In 
these systems there is no separate initiator species and a 
pair of monomers may spontaneously unite to form a living 
chain. The dynamics are thus very different and mt and Nt 
may relax on different timescales. Marques et al. stud- 
ied this class of living polymer; after linearizing dynamical 
equations they found that after a small T-jump mt relaxes 
much more rapidly than Nt which subsequently relaxes in 
time ~ N^. In numerical simulations of this same class 
of system, Milchev et al. PHI report a, 1/t decay of mean 
chain length after a time ^ N^. In the related system 
of spherical micelle aggregation |i56>.57. .58 , 59 , the num- 
ber of micelles can change and again the kinetics are very 
different to the present case. 

A summary of our findings is shown in fig. 1151 for the 
case emphasized here where a small T-jump induces a pos- 
itive velocity boost to chain growth. The figure shows the 
very different relaxation times of first and second MWD 
moments. It also shows the interesting behavior of the 
number of free initiators, ^t(O), which suffers enormous 
depletion at intermediate times despite ultimately recover- 
ing to a level very close to its initial value. We showed that 
whilst virtually all of the relaxation of the fast variables 
Nt,mt occurs on a timescale Tfast, they are actually not 
completely relaxed by this time. The final very late stages 
of their relaxation during which their values are fine-tuned 
(to relative order 1/iVoo) to the final equilibrium values. 
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Fig. 15. Time dependent recovery following a small perturba- 
tion, positive initial boost, e > ec. Top: SNt, Smt (solid line) 
become very small after Tfast while the MWD dispersion 5 At 
(dashed line) relaxes after TsIow Bottom: In the process of equi- 
libration free initiator concentration undergoes large changes 
even though its initial and equilibrium values are very close to 
one another. 



occurs on a timescale Tgiow • During this episode their val- 
ues evolve quasi-statically, enslaved to (j>t{0) according to 
6Nt9/Noo ~ —Smt/moo ~ <5(/>t(0). The novel behavior of 
^t(O), as well as its central role in late relaxation of the 
fast variables suggests this as an interesting quantity to 
measure experimentally. This might be achieved spectro- 
scopically and would have the advantage of being a rela- 
tively uninvasive probe. We remark that for bifunctional 
initiators, e.g. those used in studies by the Greer group, 
the meaning of (f>t{Q) is the concentration of "half-chains," 
i. e. those having at least one chain end which is a bare 
initiator molecule. 

We have demonstrated that the type of dynamical re- 
sponse depends on the magnitude and sign of the pertur- 
bation. For example, if the temperature change is reversed 
in sign, then the velocity boost also changes sign. Thus for 
a-methylstyrene a negative T-jump produces the positive 
velocity boost phenomenology of fig. ^| whereas a small 
temperature increase causes a negative boost. In this case, 
instead of short chains being annihilated, their number 
increases dramatically during the relaxation process. As 
far as perturbation magnitude is concerned, if the relative 

~ 1/2 

value e is less than a critical value ^ 1/Noo then the 
response is milder but remains non-linear. For a positive 
boost, the hole in the MWD is no longer complete, but 
instead has a relative depth less than unity but still much 
bigger than e. 

How do our predictions compare with experiment? In 
ref. small T-jumps were imposed on a-methylstyrene 
systems in the semi-dilute regime at concentrations high 
enough to be nearly entangled. Relaxation was monitored 
by measuring viscosity 77 as a function of time. Now gener- 
ally we expect viscosity to scale as a polymer-concentration- 
dependent power of mean chain length, rj ^ c^-ZV''', where 



the value of 7 is predicted by standard theories of polymer 
physics |6UI47j . For a general MWD, note the coefficient 
which depends on the shape of the MWD, i. e. it depends 
on the full function [^(iV)]. (This is because 77 depends not 
only on the first moment, but in general has a complex de- 
pendence on chain length distribution). Thus we predict 
an initial fast relaxation of t] lasting time Tfast — N / {9v~) 
(corresponding to the relaxation of N) followed by a much 
slower relaxation in Tgiow (this is the relaxation of the pref- 
actor C0). The fast relaxation time is independent of T- 
jump magnitude (i.e. independent of e). Now from eqs. 
(|SJ| and O we can rewrite Tfast = ^/{rk^mtot)- Using 
r — 4.7 X 10^^, mtot = Q.29gm/cm^ and a measured 
value 12 fc+ w 0.2M~^sec~^ one estimates Tfast ~ 1000 
sees. (Note that this work involved bifunctional initia- 
tors. In ref. |54| we incorrectly estimated Tfast to be twice 
the value calculated here based on a wrong estimate of 
the number of living chains which was taken to be equal 
to the number of bifunctional initiators instead of twice 
this value.) Similarly, with j2| / ~ 0.5, v~ = O.lsec"^ 
one finds TsIow ~ 1 month. Hence for these experiments 
(timescales much less than months) ij should relax in ap- 
proximately 1000 sees. In the experiments, jumps ST « 1 
"^C were imposed at various temperatures in the range 
283°K :$ T < 290°K. Now the temperature dependence 
of Tfast ~ 1/A:+ is determined by A:+ which is reported jT2] 
to vary by « 10% over this temperature range. Thus we 
expect only a slight variation in the viscosity relaxation 
time at the different temperatures studied (despite the 
fact that Noo changes significantly) . These predictions are 
very close to the experimental findings, where relaxation 
times at all temperatures were found to be about 2000 sec- 
onds, with very little variation from one temperature to 
another. This is in sharp contrast to the relative changes 
in r] itself, which varied by almost an order of magnitude. 

We have emphasized small perturbations in this work. 
In fact, the response to a large perturbation is obtained, 
qualitatively speaking, by setting e = 1 in our results. This 
case was addressed by Miyake and Stockmayer 0^ who 
considered an initial condition in which all initiators are 
free (i. e. Nq = 0). They were the first to identify the two 
timescales of the system, Tfast and TsIow, in the limit of 
very small depolymerization rates (/ <C 1). Their analyti- 
cal results applied to the t < Tfast stage where they found 
an initially uniformly translating and broadening MWD 
peak, reaching iVoo after Tfast as shown in fig. [5| Seen from 
the viewpoint of a general perturbation of magnitude e, in 
this case Tan and TsIow coincide because e = 1. That is, the 
"hole" (in this case the entire region between and iVoo) 
fills up on the same timescale in which the whole MWD 
relaxes. The linearized dynamics, which for a small pertur- 
bation appfied for t > tru due to 6(j)t{0)/4>oc,{0) becoming 
small (see section 3.3) cannot be employed. Describing dy- 
namics at times of order TsIow in this case is thus a difficult 
nonlinear problem j5 1115211^ . 

In recent experiments, Greer's group ^1^] studied 
such large perturbations using a-methylstyrene. After a 
rapid quench below the polymerization temperature, mono- 
mer concentration and MWD were analyzed after a time 
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delay. The time dependent relaxation was probed by re- 
peating the procedure for different time delays and sam- 
ples. They found different relaxation times for the mono- 
mer concentration and the MWD width, as predicted by 
theory However, Nt and rrit did not have the same 
relaxation, as mass conservation would seem to dictate. 
This may be due either to side-reactions during polymer- 
ization or living chain ionic association effects (6f..62M63| 
which have been studied theoretically Dynamic light 
scattering measurements (i 1 suggest that prior to poly- 
merization the initiators self-assemble into long polymeric 
structures and it has been suggested [12] that initiators 
may not all be equally available for polymerization follow- 
ing the T-quench. However, lifetimes of aggregate struc- 
tures would need to be extremely long for these effects to 
interfere with polymerization dynamics since as we have 
seen the MWD relaxation times are very large. 

In conclusion, our hope is that this work will motivate 
further experimental study of the dynamical sensitivity of 
living polymers to small perturbations. The drastic effect 
on the MWD in the small chain region, in particular the 
depletion or excess of free initiators, is a natural focus for 
experimental measurement. One can think of other time- 
dependent perturbations such as small amplitude thermal 
cycling which would probe interesting aspects of their ul- 
trasensitive dynamics. Finally, the systems we have ana- 
lyzed here are model starting points for the considerably 
more complex biological living polymers, actin filaments 
and microtubules, which are intimately involved in the lo- 
comotion and structural integrity of living cells |39l4l)l4l] . 
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A Derivation of Equilibrium Values moo. ^ooi 

The equilibrium MWD can be calculated by setting the 
time derivative in eq. jSJ to zero. Its solution is the Flory 
distribution, eq. where 



B Self-Consistency of Results of Subsections 
3.1 and 3.2 

B.l Coherent Chain Growth 

The solution for the Vt dynamics, eq. (|14|l . follows by solv- 
ing eq. ijnji after neglecting the (f>t{0) term on the rhs. 
For t <^ Tqs, this is a self-consistent solution since then 
the (/'t(O) term is negligible with respect to the remaining 
terms as can be seen using the solutions for vt,4>t{0) of 
eqs. inil, (UZIl, and eq. 

Eq. lfT31) is derived by substituting Vt from eq. ffT^ 
in eq. (O, replacing Dt = Doo + 5Dt, and dropping the 
5Dt term. The fact that this term can be neglected may 
be seen as follows. Even if 5Dt did not decrease but re- 
mained of order 5Dq throughout this regime this would 
lead to replacing D^o by I?oo(l + e) in all results obtained 
in section 3 (see eq. CSJ)- This upper bound on the ef- 
fect of 5Dt can easily be seen to lead to higher order in 
e terms (see eg. eq. (jlTIl ') and the 6Dt term can therefore 
be neglected. 

Now let us consider the propagator of the dynamics 
of eq. (|f 5|l . which appears in eq. H16|l . For short times, 
t ^ Tfast, one may replace Vt in eq. (|15|l by vq. Eq. (|15|l 
then has constant coefficients and its propagator including 
the boundary condition may be calculated exactly (HH] : 

I 



^^-(N-N' -votf /(AD^t) _^ ^-N'vo/D^^~(N+N'-votf/{4:D^t) 

Here erfc(a;) = 1 — erf (a;), where erf is the error function. 
The solution for 0t(O) is derived by substituting eq. 
in eq. lfTC|l and setting iV = 0. One has 



■ erfc(a;) + 2x^eYic{x) ~ 2TT-^''^xe-^^ 



1/2 



(i<Tfast) (43) 



Solving the system of eqs. 
has 



1 + / + 



(HOI) and 



© for 



(40) 



, iVoo one 




(l-/)2 + r(I + 3/) + 



1/2 > 



"^tot 



Here we used the fact that since the values of N' con- 
tributing to the integration are much smaller than iVoo, 
one may approximate 4>o{N') by (/)oo(0) in the integrand 
of eq. (|f 6(1 (with relative errors of order 1/iVoo). Eq. 1(17(1 of 
the main text follows from eq. 1(43(1 by taking the two cor- 
responding limits. Note that eq. ((43(1 approximates (/'t(O) 
for all t <C Tfast, including times of order t* . 



(41) 



For r ^ 1 and assuming (1 — f )/r^^^ ^ 1, by expanding 
eq. ((41(1 in powers of r and keeping only the first term of 
the expansion one recovers eq. ^ . The values of Woo , D^o 
in eq. ^ are obtained by substitution of moo , -^oo in eq. 



B.2 Hole Filling 

Eq. ((18(1 is a self-consistent solution of eq. ^ since for 
t 3> Tqs, assuming the validity of eqs. ((18(1 and ((17(1 one 
may verify that eq. jnj is satisfied with the term on the 
Ihs being much smaller in magnitude than both terms on 
the rhs which balance one another. 
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Consider now the validity of eq. ((11^ for Tfast ^ t <C 
Tfiii. Substituting eq. (|21|l into eq. H18(l one has explicit 
solutions for rrit, Nt,Vt, Dt- One thus sees that in the time 
regime of consideration the velocity term in eq. © is ex- 
ponentially small and therefore can be deleted in eq. H19|l 
with very small error. 

The MWD dynamics are therefore described by eq. 
H19|) whose propagator satisfying reflecting boundary con- 
ditions is 



Gf"(7V,iV') = 



1 



|g-(W-W')V(4Doot) _,_ g-(JV+JV')V(4D„„t)| , (44) 

The solution of eq. 119|l is thus 



diV>*„_(7V')G^l_(iV,iV') 



dN' MN' - eN^/e)Gf'''{N,N') 

(Tfast < t < Tfiii) . (45) 

where Across Satisfies Tfast < Across < Tqs, i. 6. it belongs 
to the time regime described simultaneously by eqs. H15|l 
and (|19|l . In going from the first integral in eq. l|l5l) to 
the second we used the fact that at tcross the MWD has 
approximately the shape of the initial MWD, boosted in 
the positive direction by eNao/0. The effect of diffusion 
during the initial boost of the MWD may be shown to be 
small. Notice also that evolution in the last expression in 
eq. (|45|l starts at t = 0; for t ^ icross we may to leading 
order replace t — Across in C*'^ in eq. by t. 

Substituting eq. (|44|l in eq. (|45() and setting = one 
obtains 



«(0) 



rfc 



Cf'f 



(Tfast < t < Tfill) (46) 



Considering times much shorter than the crossover time 
Tfill, eq. is recovered. 



C Linearized Dynamics, t ^ Tfiu 

Ultrasensitivity is such that the dynamics cannot be lin- 
earized until the later stages, t 3> Tfin. For these times we 
may linearize eq. lO by dropping terms proportional to 
products of (5(/>t, 6vt, SDf. One has 



dt 



ON 



6vt 



(47) 



with reflecting boundary conditions at the origin. In eq. 
(|47() the 5Dt term was neglected since its magnitude is 
smaller by l/N^o than the corresponding Svt term. The 
(5-function in the source term in eq. H47|) arises from lin- 
earizing the boundary condition in eq. ij^l. Notice that 



by definition, (5(/'t(0) is normalized to zero and that the 
source term in eq. (|47ll preserves this normalization since 
its integral over all N is zero. 

Performing the corresponding linearization in eq. © 
one has 



-rSvt 
dt 



1 



1 



Svt DooSMO) 



Tfast 



Tfast 



(48) 



Now since eq. (|47|l is of the same form as eq. (|15|) . 

its propagator, G\^^°^'^{N, N'), is given by eq. ig^l with 
Vq replaced by v^o- Thus if (5(/>tj^ (N) is known at time i^,, 
then the solution of eq. (|47|l for subsequent times is 



SMN) « / dN'G'rtT{N,N')5c^tAN') + Rt{N) , 
Jo 

(49) 

where 



Rt{N)= dN' dt'G';^1f'{N,N')fit'{N') 



' dt' [0^(jV)-Gl'— (A^,0)] • (50) 



Here the last expression for Rt (N) is obtained by perform- 
ing the N' integration using the expression for /if from eq. 
(I47|l and the fact that (poo{N) remains unchanged when 
evolved with G""'=^''. In eq. ^ we also used 



Svt ~ -D^dMO) 



(51) 



which is valid throughout the linearized time regime as 
one may derive using eq. H48|) . The solution of eqs. (|47|l 
and (|48|l is thus found by setting = in eq. 1)49(1 , solving 
for SiptiO), and then using the calculated expression in eqs. 
(gni and (EU to find <5<^t(A^) and Svt- Let us now perform 
this analysis for two time regimes. 

(i) Tfill ^ i ^ Tgiow In this regime, we will see the 
source term can in effect be ignored. The top expression 
on the rhs of eq. (|22() is a self-consistent solution of eq. (|49|l 
which is proved as follows. As one may see from eqs. (|21|l 
and (|22|) there exists a constant f3 of order unity such that 
for t > /Sthii the relative perturbation in iptiN) is much 
smaller than unity and eq. l|T7)l applies. (More precisely, 
for any desired relative smallness there exists a different 
constant /3.) Thus setting ti, = (3Tf^\\ and A^ = in eq. 
and using the expression for 5<t)t (0) from eq. 1(22(1 one finds 
that with appropriate choice of J, eq. 1(49(1 is satisfied with 
Rt (0) being much smaller than the remaining two terms in 
eq. ((49(1 . The integral term in eq. 1(49(1 is calculated using 
the fact that as seen from eqs. ((21II . (I22II . and fig. 
5(t)tdN) is of order (5(/>tjO) > e<Poo{N) for A^ < eN^/O 
and of order e(j)oo{N) for larger A^. We also used 



Gn'"(0,iV') 



l/[D^{t-t')f' {N' ^ D^{t - t')) 
{N' :^ Doo{t - t')) 

(t,t'<Tslow) (52) 



which may be derived from eq. ((42(1 . 
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(ii) t » Tsiow Using eq. (|^ one finds 



linear 



G 



where 



(0,iV') «5oo + (5t-5oo)A(7V') (i»r,iow), (53) 



St = G'r^'iO, 0) , Soo= cl^ociO) = l/^oo , 



X{N') 



N' 
2nZ 



(54) 



Now setting II — TsIow in eq. I^?^ and using eqs. 
and one has 



bt{Q) ^{St-S^) / dN' \{N')5^r^,^^{N') + 

JQ 



1 



45', 



oo'^slow JTsI, 



dt' {St-f -S^)6cPt,{0) (55) 



where we used the fact that S(I)t^-^^^{N') is normahzed to 
zero. 

Now the self-consistency of the t > Tgiow expression on 
the rhs of eq. (|22(l is proved by substituting eq. (|22(l in eq. 
H55|) . Since Scpr^ia^iN') is of order e/Noc, the magnitude 
of the integral of S(j)r^^^^ {^') in eq. 1)55(1 is of order e. The 
last integral term in eq. ((55|l is evaluated using eq. H22|l 
and 

St >S'oo ~ 

165oo7r-i/2 (T,iow/i)'/' e-*/(i6-=.ow) » ^^j^^) 

(56) 

which may be proved using eqs. (|^ and One finds 
that for t 3> Tgiow all terms in eq. (|55|) have the same 
time dependence and are of the same order of magnitude. 
Thus eq. 1(55(1 may be satisfied by appropriate choice of 
the numerical coefficient M in eq. ((22(1 . 



D Self-Consistency of the Results of Section 
5 

D.l Coherent Chain Shrinking 

The validity of eq. ((14(1 is verified using the same argu- 
ments as in the first paragraph of Appendix B but now 
using the expression for (j>t (0) from eq. I(28() . 

Setting iV = in eq. one obtains eq. (|^ with x 
now replaced by — x. Considering times much greater and 
much less than t* leads to eq. ((28(1 . 



D.2 Peak Decay 

The validity of eq. I(18() is proved using exactly the same 
arguments as the ones in the first paragraph of the corre- 
sponding section of Appendix B, but now using the expres- 
sion for (ptiO) from eq. ((291) . One finds that the crossover 



time Tqs ~ Tfast lii(iVoo) is different from the respective 
crossover time found in Appendix B. 
The solution of eq. 1(19(1 is now 



MN)- / dN'MN' -eN^/0)Gf*'{N,N') + 
Jo 

Gf"" {N, 0) / dN' 00 {N') (Tfast « t « r,iow) (57) 

Jo 

Here the evolved MWD consists of two parts, correspond- 
ing to the two terms on the rhs (see fig. Illf c)). The first 
is the initial MWD, (po, shifted by eNoo/d towards smaller 
N . The second part is the excess peak accumulated by 
the negative boost at the origin and whose amount was 
the polymer initially lying between zero and eNoo/0. Here, 
as in section 3, the effects of diffusive broadening during 
t < Tfast have been neglected, and Tfast in eq. 1(57(1 has been 
replaced by zero. 

Now substituting eq. 1(44(1 in eq. ((57(1 and setting N — 
one obtains eq. ((29() similarly to the derivation of eq. ((21(1 
in Appendix B. 



E Self-Consistency of the Results of Section 6 

For times t <^ Tfast, all results and analysis are identi- 
cal to those for the strong perturbation (e > Cc) case in 
the regime t ^ t* . Eq. 1(32(1 can be shown to be a self- 
consistent solution using the same arguments as those of 
Appendices B and D to prove eqs. ifT^ and ifT^ . respec- 
tively. For times longer than Tfast hncarization of the dy- 
namics can be performed and the results of Appendix C 
directly apply. 
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